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Abstract 

We describe a wide class of polynomials, which is a natural generalization of Hurwitz stable poly- 
nomials. We also give a detailed account of so-called self-interlacing polynomials, which are dual to 
Hurwitz stable polynomials but have only real and simple zeroes. All proofs are given using proper- 
ties of rational functions mapping the upper half-plane of the complex plane to the lower half-plane. 
Matrices with self-interlacing spectra and other applications of generalized Hurwitz polynomials are 



discussed. 
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Introduction 

The polynomials with zeroes in the open left half-plane of the complex plane are called Hurwitz stable, 
or just stable. They play an important role in several areas of mathematics and engineering such as stability 
and bifurcation theory, control theory etc. The Hurwitz stable polynomials are well-studied and closely 
related to the moment problem, spectral theory, operator theory and orthogonal polynomials. The Hurwitz 
stable polynomials were studied by C.Hermite, E. Routh, A. Stodola, A. Hurwitz [H] [35l [I7J [22 [12] . The 
present paper was motivated by some problems of the bifurcation theory and concerns the most natural 
(from our point of view) generalization of Hurwitz stable polynomials in the following way. 

Let 

p( z ) =Po(z 2 ) + zp 1 {z 2 ) 

be a real polynomial, where po and p\ are its even and odd parts, respectively. Consider its associated 
rational function 

Po(u) 

In the present paper we study polynomials whose associated function $ is a so-called i?-function, that is, 
the function mapping the upper half-plane of the complex plane to the lower half-plane. These functions 
have real, simple, and interlacing zeroes and poles [HI [16]. Such functions play a very important role 
in spectral theory of self-adjoint operators [37] , in the moment problem [3] |TJ [S3] , in the theory of the 
Stieltjes string [13], in the theory of oscillation operators 10 , in the theory of Jacobi and Stieltjes continued 
fractions [41] etc. In particular, it is very well known [22j [3 [l 2 . 6 that the polynomial p is Hurwitz stable if 
and only if its associated function $ is an R- function with negative zeroes and poles, and degp ^ degpi. 
In Section[3jof the present paper, we prove this fact using methods of complex analysis, and then from this 
fact we obtain all main properties of Hurwitz stable polynomials and the main criteria of the polynomial 
stability including connections with total positivity and Stieltjes continued fractions. Section 1X21 is devoted 
to polynomials with zeroes in the closed half-plane of the complex plane. Note that in [3, 20 there was 
proved that the infinite Hurwitz matrix of a quasi-stable polynomial is totally nonnegative. In Section [3. 2i 
we proved the converse statement, which is new. 

In Section [4] we investigate polynomials whose associated function $ maps the upper half-plane to 
the lower half-plane (with degpo ^ degpi) but has only positive zeroes and poles in contrast with the 
Hurwitz stable case . We call those polynomials self-interlacing. If p(z) is a self-interlacing polynomial, 
then it has only simple real zeroes which interlace with zeroes of the polynomial p(—z). More exactly, if 
Xi, i — 1, . . . , n, are the zeroes of the self-interlacing polynomial p (of type I), then 



< Ai < -A 2 < A 3 < . . . < (-1)"-%. 



This fact is very curious, because the deep structure of self-interlacing polynomials is very close to the 
structure of the Hurwitz polynomials but the former ones have only real and simple zeroes. The main 
result that reveals this connection is Theorem 14.81 which says that a polynomial p(z) = po(z 2 ) + zpi(z 2 ) 
is self-interlacing if and only if the polynomial q(z) = p(— z 2 ) — zpi(z 2 ) is Hurwitz stable. Section [4] is 
devoted to the comprehensive description of self-interlacing polynomials and includes analogues of the 
Hurwitz criterion of stability and of the Lienard and Chipart criterion of stability, and some other criteria 
of self-interlacing. We also give some examples of matrices whose characteristic polynomials are self- 
interlacing and describe some numerical observations related to the self-interlacing polynomials. Note that 

there are three i?-functions connected to the self-interlacing polynomials. Namely, let p be self-interlacing. 

p' 
Then it has only real and simple zeroes, so its logarithmic derivative — is an i?-function (see e.g. |16j). 

Moreover, p is self-interlacing if and only if its associated function $ = — is an i?-function with positive 

Po 

zeroes and poles (Theorem I4.3[) . Finally, p is self-interlacing if and only if the function is an 

p(-z) 



R- function (see the proof of Theorem I4.3|) . 

The full generalization of the Hurwitz stable polynomials (in terms of the function $) is given in Sec- 
tion [5] We describe the polynomials whose associated function $ is an R- function without any restriction 
on the signs of its zeroes and poles. We call such polynomials generalized Hurwitz polynomials. More 



exactly, the polynomial p of degree n is called generalized Hurwitz of order k (of type I) , $J k $J 
if it has exactly k zeroes in the closed right half-plane, all of which are nonnegative and simple: 

s^ fii < ^ 2 < • •• < Hk, 

such that p(—/m) t^ 0, i = 1, . . . , k, and p has an odd number of zeroes, counting multiplicities, on each 
interval (— /Ufc, — Atfc-i), • • • > (— A*3, —^2), (~A*2, — Hi)- Moreover, the number of zeroes of p on the inter- 
val (— /ii,0) (if any) is even, counting multiplicities. The other real zeroes lie on the interval (—00, — //&): 
an odd number of zeroes, counting multiplicities, when n — 21, and an even number of zeroes, counting 
multiplicities, when n = 21 + 1. All nonreal zeroes of p (if any) are located in the open left half-plane of 
the complex plane. 

We prove a generalization of Hurwitz theorem, Theorem 15.61 which, in fact, is an analogue of Hurwitz 
criterion of stability for the generalized Hurwitz polynomials. We also establish a generalization of the 
Lienard and Chipart criterion of stability, Theorem l5.8l Note that in [21] G. F. Korsakov made an attempt 
to prove this generalization of the Lienard and Chipart criterion. However, his methods did not allow him 
to prove the whole theorem (for details, see Section IS - !"]) . 

We also should note that V. Pivovarchik in [31] describes some aspects of the generalized Hurwitz poly- 
nomials, however, in different form, replacing the lower half-plane by the upper half-plane. In particular, 
he indirectly establishes that the polynomial p is generalized Hurwitz (shifted Hermite-Biehler with sym- 
metry, in Pivovarchik's terminology) if and only if its associated function $ is an i?-function. Precisely, 
the cases considered in [31] are degpo > degpi and degpo = degpi with the leading coefficients of po and 
pi of the same sign. Neither generalizations of Hurwitz theorem nor generalizations of the Lienard and 
Chipart criterion, nor the theory of self-interlacing polynomials, nor connections of the generalized Hurwitz 
polynomials with continued fractions were established in |31j . However, V. Pivovarchik with H. Woracek 
described [33] [34] a generalization of his shifted Hermite-Biehler polynomials with symmetry (which are 
generalized Hurwitz polynomials up to rotation of the upper half-plane to the left half-plane) to entire 
functions. The main result of [33; 34 is an analogue of our Theorem 15.31 for entire functions, provided by 
two different methods. The method in [34] coincides with that of [31], where as [34] uses the theory of de 
Branges and Pontryagin spaces. Our proof of Theorem [231 differs from the methods used by V. Pivovarchik 
and H. Woracek and is based on using the properties of ^-functions. 

It is curious that V. Pivovarchik came to his shifted Hermite-Biehler polynomials and entire functions 
with symmetry due to his investigation of spectra of quadratic operator pencils [30 and Sturm-Liouville 
operators whose boundary conditions contain spectral parameter |32| I40j . The last two articles deal with a 
generalization of the Regge problem, which has a direct relation to the scattering theory. So the distribution 
of zeroes of the generalized Hurwitz polynomials has a physical interpretation. 

Finally, in Section [5.21 we explain how the generalized Hurwitz polynomials can be used in bifurcation 
theory. Also, we describe how specific zeroes of a real polynomial p(z, a) move with the movement of the 
parameter a along R if p(z, a) is generalized Hurwitz for all real a. 

1 Auxiliary definitions and theorems 

Consider a rational function 

R(z) = «$-, (1.1) 

where p and q are polynomials with complex coefficients 

p(z) = a z n + aiz 71 ^ 1 + ■ ■ ■ + a n , a , 01, . . . , a n £ R, a > 0, (1.2) 

q(z) = b Q z n + biz"- 1 + ... + b n , b ,bi,...,b n eR, (1.3) 

so degp = n and degq ^ n. If the greatest common divisor of p and q has degree I, < I < n, then the 
rational function R has exactly r = n — I poles. 



1.1 Hankel matrices and Hankel minors 

Expand the function (jl.ip into its Laurent series at oo: 



R(z) 



Here s ? - = for j < n — 1 — m and s„_!_ 



so 
S-] H h — + — 



£l 

z 2 



£2 

Z 3 



(1.4) 



&o 
a ' 



where m = deg q. 



The sequence of coefficients of negative powers of z 

s ,si,s 2 , ■ • ■ 
defines the infinite Hankel matrix S:=S(R) := \ Si+j||°° = o ■ 
Definition 1.1. For a given infinite sequence (sj)J^ , consider the determinants 



Dj(S) 



so 

•si 

S.7-1 



Si s 2 
S2 S 3 



Sj-i 



S2J-2 



.7 = 1,2,3, 



(1.5) 



i.e., the leading principal minors of the infinite Hankel matrix S. These determinants are referred to as 
the Hankel minors or Hankel determinants. 

An infinite matrix is said to have finite rank r if all its minors of order greater than r equal zero 
whereas there exists at least one nonzero minor of order r. Kronecker 24 proved that, for any infinite 
Hankel matrix, any minor of order r where r is the rank of the matrix, is a multiple of its leading principal 
minor of order r. This implies the following result. 

Theorem 1.2 (Kronecker [24] ) . An infinite Hankel matrix S = \\si+j\\^ has finite rank r if and only if 



D r (S) ^ 0, 

DjiS) = 0, for all j > r. 



(1.6) 
(1.7) 



The following theorem was established by Gantmacher in [T2] , 
Theorem 1.3. An infinite Hankel matrix S — ||s-i+j||o° has finite rank if and only if the sum of the series 

so 



R{z) 



£i 

..2 



£2 



is a rational function of z. In this case the rank of the matrix S is equal to the number of poles of 
the function R. 



Theorems 11.21 and 11.31 imply the following: if the greatest common divisor of the polynomials p and q 
defined in (|1.2[) - (|1.3|) has degree I, ^ I ^ to, then the formulae (|1.6|) - (|1.7|) hold forr = n-l for the 
rational function (|l.l|l . 

Let Dj(S) denote the following determinants 



D 3 (S) 



Sl 


S2 


S3 


S'2 


S3 


S'4 


■ s .; 


Sj+1 


Sj+2 



Sj+1 



S2J-1 



j = 1,2,3,.. 



(1.8) 



With a slight abuse of notation, we will also write Dj(R) := Dj(S(R)) and Dj{R) :=Dj(S(R)) if the matrix 
S = S(R) is made of the coefficients (jl.41) of the function R. 

Theorems II .21 and II. 31 have the following simple corollaries, which will be useful later. 



Corollary 1.4 (|16j). A rational function R represented by the series (jl.4p has at most r poles if and only if 

Dj(R) = for j = r + l,r + 2, . . . . 

Corollary 1.5 f|16j). A rational function R with exactly r poles represented by the series (|1.4p has a pole 
at the point if and only if 

D r -i(R) ^ and D r (R) = 0. 

Otherwise, D r (R) ^ 0. 

In the sequel, we are also interested in the number of sign changes in the sequences of Hankel minors 
Dj{R) and Dj(R), which we denote b£| V(D (R),Di(R), ... ,D r (R)) andV(D (R),Di(R),...,D r (R)), 
respectively. These numbers exist, since R is a real function. 

In his remarkable work [S], G. Frobenius proved the following fact that allow us to calculate the number 
of sign changes in the sequence Do(R),Di(R), . . . , D r {R) in the case when some minors Dj(R) are zero. 

Theorem 1.6 (Frobenius [SHU]). If, for some integers i and j (0 < i < j), 

Di(R)?0, D i+l {R) = D i+2 (R) = ■ ■ • = Di+j(R) = 0, D i+j+l {R) £ 0, (1.9) 

then the number Y(D (R), Di{R) 1 D2(R), ■ ■ ■ , D r (R)) of Frobenius sign changes should be calculated by 
assigning signs as follows: 

signA+,(i?) = (-l) iii ^signA(i?), v = l,2,...,j. (1.10) 

Definition 1.7. A matrix (finite or infinite) is called totally nonnegative (strictly totally positive) if all 
its minors are nonnegative (positive). 

Definition 1.8. An infinite matrix of finite rank r is called m-totally nonnegative (m-strictly totally 
positive) of rank r, m $J r, if all its minors are nonnegative (positive) up to order m inclusively. 

An r-totally nonnegative (r-strictly totally positive) matrix of rank r is called totally nonnegative 
(strictly totally positive) of rank r. 

Definition 1.9. An infinite matrix is called positive (nonnegative) definite of rank r if all its leading 
principal minors are positive (nonnegative) up to order r inclusive. 

For the matrix A (finite or infinite), we denote its minor of order j constructed with rows ii, i 2 , ■ ■ ■ ,ij 
and columns lx, l 2 , . • . , lj by 

A (? f ? 

\ii H ... lj 

Definition 1.10 ([10 J. An infinite matrix A of finite rank r is called m-sign regular, m ^ r, if all its 
minors up to order m (inclusively) satisfy the following inequalities: 

Yl l k+ J2 l k , fii i 2 ... ij 



(-i)& *=i aij £ ;;; gi>o, j = i,...,m. (1.11) 

An r-sign regular matrix of rank r is called sign regular of rank r. 

In |12j there was proved the following criterion of total positivity of infinite Hankel matrices of finite 
order. 

Theorem 1.11. An infinite Hankel matrix S = ||si+j'||o° of finite rank r is strictly totally positive if and 
only if the following inequalities hold: 

J j = l,...,r. (1.12) 



1 Here we set D (R) = D (R) = 1. 



In Pj5], the following fact was proved. 

Theorem 1.12. An infinite Hankel matrix S = ||sj+j||g° of finite rank r is sign regular if and only if 
the following inequalities hold: 



Dj(S)>0, 

(-im-(s) >o, 

1.2 Matrices and minors of Hurwitz type 



3 = 1, 



(1.13) 



In this section we introduce infinite matrices (matrices of Hurwitz type) associated with the polynomials 
p and q defined in (ll.2[) - (ll.3[) and discuss the Hurwitz formulaaethat connect Hankel matrices with matrices 
of Hurwitz type. 

Definition 1.13. Given polynomials p and q from ()1.2j) 
if degg < degp, that is, if b = 0, then 



H(p,q) 



if degg = degp, that is, bo ^ 0, then 



H(p,q) = 



(ao 


a\ 


0-2 


03 


«4 


05 . . . 





In 


b-2 


h 


bi 


&5 ... 





'1-0 


ax 


a-2 


f'3 


04 








bi 


b 2 


h 


& 4 ... 


V : 












(bo 


h 


b 2 


b,, 


bi 


h ..- 





«o 


a i 


o 2 


0:>, 


a4 





bo 


hi 


b 2 


h 


o 4 ... 








no 


a i 


a 2 


a 3 ... 



1.3[) . define the infinite matrix H (p, q) as follows: 

\ 

(1.14) 

) 
\ 



(1.15) 



The matrix H(j>, q) is called an infinite matrix of Hurwitz type. We denote the leading principal minor 
of H(p, q) of order j, j = 1, 2, . . ., by rjj(p, q). 

Remark 1.14. The matrix H(p, q) is of infinite rank since its submatrix obtained by deleting the even or 
odd rows of the original matrix is a triangular infinite matrix with ao ^ on the main diagonal. 
Together with the infinite matrix H(p,q), we consider its specific finite submatrices: 

Definition 1.15. Let the polynomials p and q are given by (|1.2j) - (|1.3[) . 
If degg < degp = n, then we construct the In x 2n matrix 



/ 


'bi b 2 b 3 ... b n 


... 





°\ 




o-o a-i 0,2 ... a n -i a n 


... 










bi b 2 ... o n _i b n 


... 








H2n(p,q) = 


ao ai ... a n _2 fln-i 


a„ ... 










0... ao ai 


a 2 ... 


a„ 







0... bi 


b 2 ... 


b„ 





\ 


K ... a 


ai ... 


a n -i 


a n J 


n, then we construct the (2n+l) x (2n+l) matrix 








(a ai a 2 ... a n -x a„ 


.. 





°\ 




b h b 2 ... b n -i b n 


.. 










ao ai ... a n _2 in-i 


a n . . 








i 2n +i(p,q) = 


b h ... &„_ 2 &«-i 


b n ■■ 










0... ao ai 


a 2 ■■ 


a„ 







... b bi 


b 2 .. 


b n 







\0 ... a 


ai .. 


dn-i 


a n J 



(1.16) 



(1.17) 



Both kinds of matrices 'H.2n(Pi o) an d %2n+i(Pj q) are called finite matrices of Hurwitz type. The leading 
principal minors of these matrices will be denoted bj|3 V^(p, g). 

In his celebrated work [IT], A. Hurwitz found relationships between the minors Dj(R) and Dj(R) 
defined in (|1.5p and (|1.8[) and the leading principal minors ??;(j>, q) of the matrix H(p, g) (see Dcfinition ll.l3|) . 

Lemma 1.16 ( [T7J [2H HH [SI US] ) ■ Let the polynomials p and q be defined in (|1.2l) - (|1.3p and let 



q(z) so si 

— - = S-l + — + -j 

p[z) z z 



The following relations hold between the determinants t]j (p, q) and Dj(R),Dj{R) defined in (|1.5[) and (jl.81) , 
respectively. 



Ifdegq < degp, then 



If degq = degp, then 



r,2 j (p,q) = a 2 j D j (R), j = l,2,...; (1.18) 

m J+ i(p,q)^(-^y^ +1 Dj(R), j = o,i,2,... (i.i9) 



7 l2j+1 ( Pl q)^b a 2 'D J (R), j =0,1,2,...; (1.20) 

7 72j (p, ( j) = (-l)^ 1 6 oa ^- 1 5 J _ 1 (i?), j = 1,2,... (1.21) 

From (fLTH (fTTTj) we have 

If degg < degp, then 

V l (p, 9 ) = ao 1 7 ?J+1 (p,g), i=l,2,...,2n; (1.22) 

If degg = degp, then 

V l (p,g) = 6 ( 7 1 7 ?!+1 (p, (Z ), i=l,2,...,2n + l. (1.23) 

1.3 Rational functions mapping the upper half-plane to the lower half-plane 

All main theorems related to the generalized Hurwitz polynomials are based on properties of the 
following very important class of rational functions: 

Definition 1.17. A rational function F is called i?-function if it maps the upper half-plane of the complex 
plane to the lower half-planqf|: 

lmz>0^ ImF(z) < 0. 

By now, these functions, as well as their meromorphic analogues, have been considered by many authors 
and have acquired various names. For instance, these functions are called strongly real functions in the 
monograph [36 due to their property to take real values only for real values of the argument (more general 
and detailed consideration can be found in [7], see also [IS])- 

The following theorem provides the most important (for the present study) properties of i?-functions. 
Some parts of this theorem was considered in j29] [22] [7j [T21 H] |6] . An extended version of this theorem can 
be found (with proof) in [HI Theorem 3.4]. 

Theorem 1.18. Let h and g be real polynomials such that deg/i — 1 $J degg = n. For the real rational 
function 

F= h - 
9 

with exactly r (^ ri) poles, the following conditions are equivalent: 

1) F is an R-function: 

Imz> O^ImF(z) < 0; (1.24) 



2 That is, Vj(p, q) is the leading principal minor of the matrix "H2n(p,q) of order j if degg < degp. Otherwise (when 
■gq = degp), Vj(p, q) denotes the leading principal minor of the matrix "H2n+l(p, q) of order j. 
3 In |16l these functions are called R-functions of negative type. 



2) The function F can be represented in the form 



where 



F(z) = -az + [5 + V — — , a)OJ,WjeM, (1.25) 

^— : f z + u)j 

7 J - = ^r4>0, i = l,...,r; (1.26) 



3) T/ie zeroes of the polynomials g(z) and h(z), where g = 3//, ft- = h/f, f = gcd(g, ft), are rea/, simple 
and interlacing, that is, between any two consecutive zeroes of one of polynomials there is exactly one 
zero of the other polynomial, counting multiplicity, and 

3oj eR : g(u)h'{uj) ~ g'(u>)h(cu) < 0; (1.27) 

4) The polynomial 

g{z) = Xp(z) + m {z) (1.28) 

has only real zeroes for any real A and \x, A + \r 7^ ; and the condition (|1.27[) is satisfied; 

5) Let the function F be represented by the series 

F(z) = B_ 2Z + 5_! + - + ^ + ^ + ..., (1.29) 

z z z Z A 

then s_2 ^J and s_i € M. TTie following inequalities hold 

D J (F)>0, i = l,2,...,r, (1.30) 

where determinants Dj(F) are defined in (|1.5|) . 

Remark 1.19. If a function i 7, maps the upper half-plane to the upper half-planq^|, then the function — _F 
is an i?-function. 

Note that from (|1.25p it follows that if F is an i?-function, then the function F is decreasing between 
its poles, so — F is obviously increasing. Therefore, from Theorem ll.lSl we obtain the following well-known 
fact. 

Corollary 1.20. Let g and h be real polynomials such that \ degh — deg<?| ^ 1 and let the zeroes of the 
polynomials g and h be real and simple. If zeroes of g and h interlace, then the function — is decreasing 

g 

or increasing between its poles. 

Remark 1.21. If R defined in (11.11) is an i?-function, then by Theorem 11.181 we have degq ^ degp — 1, 
that is, the equality bo = implies b\ ^ or, in other words, 6q + b\ ^ 0. 

In the sequel, we also use the following fact [16], which is a simple consequence of a theorem of 
V.Markov [7]. 

Theorem 1.22. Let p and q be real coprime polynomials such that \ degp — degg| ^ 1. If the function 
R = q/p is an R-function, then the functions Rj = g") /pti) ; j = l j . , . ) degp — 1, are also R-functions. 



The following simple properties of R- functions can be established using Theorem II. 181 
Theorem 1.23. Let F and G be R-functions. Then 

1) the function R = F + G is an R-function; 

2) the function —F(—z) is an R-function; 

3) the function — is an R-function; 



'In 16 such functions are called .R-functions of positive type. 



4) the function zF(z) is an R-function if all poles of F are positive and F(z) — > as z — > oo. 

Now consider again the function R denned in (|1.1|) , By (jl.2|) - (jl.3|l . the degree of its numerator is not 
greater than the degree of its denominator. If R is an i2-function, then one can find the numbers of its 
negative and positive poles (see e.g. [16]). 

Theorem 1.24. Let a rational function R with exactly r poles be an R-function of negative type and let 
R have a series expansion (|1.29[) . Then the number r_ of negative poles of R equal$j 

r- =V(l,D 1 (R),D 2 (R),...,D k (R)), (1.31) 

where k = r — 1, if R(0) = oo, and k = r, if |i?(0)| < oo. The determinants Dj(R) are defined in (TIT 



As was shown in Q2] (see also [11]), if R is an i?-function has exactly r poles, all of which are of the 
same sign, then all minors Dj(R) are non-zero up to order r. 

Corollary 1.25 ([12j). Let a rational function R with exactly r poles (counting multiplicities) be an 
R-function of negative type. All poles of R are positive if and only if 

D j (R)>0, j = l,...,r, 

where the determinants Dj(R) are defined in (jl.8p . 

Corollary 1.26 (|16j). Let a rational function R with exactly r poles be an R-function of negative type. 
All poles of R are negative if and only if 

(-l) j D j (R)>0, j = l,...,r, 

where the determinants Dj(R) are defined in (|1.8[) . 

One can use (|1.18p - (|1.23p to obtain criteria for the function R defined in (jl.ip to be an iZ-function in 
terms of Hurwitz minors. 

Definition 1.27. For a sequence of real numbers ao, oi, . . . , a„, we denote the number of sign changes in 
the sequence of its nonzero entries by v(ao, oi, . . . , a„). The number v(ao,ai, . . . , a n ) is usually called the 
number of strong sign changes of the sequence ao, Oi, • • • , a n - 

In [IB], the following theorem was established. 

Theorem 1.28. Let R be a real rational function as in (|1.1|) . If R is an R-function with exactly n poles, 
that is, gcd(p,q) = 1, then the number of its positive poles equals v{ao,a\, . . . , a n ). In particular, R has 
only negative poles if and only i/|£| a,j > for j = 1, 2, . . . , n, and it has only positive poles if and only if 
aj-iaj < for j = 1, 2, . . . , n. 

Remark 1.29. It is easy to see that Theorem 1 1.281 is also true for an iZ-function with degree of numerator 
greater than degree of denominator. 

For R- functions with only negative poles, there are a few more criteria [16] . 

Theorem 1.30. The function (jl.ip . where dcgq < dcgp, is an R-function and has exactly n negative 
poles if and only if one of the following equivalent conditions holds 

1) a n > 0,a„_i >0, ...,a > 0, Vi(p,g) > 0,V 3 (p,q) > 0, . . . , V 2n -i (p, q) > 0: 

2) a n >0,b n >0,6 n _x >0, ...,6i >0, Vi(p,q) >0,V 3 {p,q) > 0, . . . , V 2 „-i(p,q) > 0; 

3) a„ > 0, a n _i > 0, . . . , a > 0, V 2 (p, q) > 0, V 4 (p, q) > 0, . . . , V 2 „(p, q) > 0; 

4) a„ > 0,6„ > 0,6 n _i >0, ...,6i > 0, V 2 (p,q) > 0,V 4 (p,g) > 0, . . . , V 2 „(p,g) > 0; 



5 Recall that the number V(l, Di(R), L>2(R), . . . , D^{R)) of Frobcnius sign changes must be calculated according to Frobe- 
nius rule provided by Theorem ll.6l 

6 In fact, the coefficients must be of the same signs, but we assume that ao > (see HI. 21 0. 
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where V$(p, q) are defined in Definition ] 1.151 

In the case degg = degp we have similar criteria. 

Theorem 1.31. The function (|1.1[) . where degq(z) — degp(z), is an R-function and has exactly n negative 
poles if and only if one of the following equivalent conditions holds 

1) a n > 0,a n _i >0, ...,a >0, V 2 (p,q) > 0,V 4 (p,g) > 0, . . . , V 2n (p, q) > 0; 

2) a„ > 0, 6 n > 0, 6„_i >0,...,6 >0, V 2 (p,«) >0,V 4 (p,g) > 0, . . . , V 2n (p,g) > 0; 

3) a„ > 0,a n _i >0, ...,a > 0, Vi(p,g) > 0,V 3 (p,g) > 0,.. . , V 2n+ i(p,q) > 0; 

4) a„ > 0,6„ > 0,6 n _i >0,...,6 > 0, Vi(p,g) > 0,V 3 (p,g) > 0, . . . , V 2n +i (p, q) > 0; 

where Vj(p, g) are defined in Definition ] 1.151 

At last, we recall a relation between ^-functions with nonpositive poles and some properties of Hurwitz 
matrices |16j . 

Theorem 1.32 (Total Nonnegativity of the Hurwitz Matrix). The following are equivalent: 

1) The polynomials p and q defined in (|1.2|) - (|1.3p Ziaue on/y nonpositive zeroes (or q(z) = 0), and the 
function R defined in (|1.1[) is either an R-function (or R(z) = 0). 

2) The infinite matrix of Hurwitz type H(p 1 q) defined in (|1.14[) - (|1.15[) is totally nonnegative. 

Theorem ll. 321 implies the following corollary |16) . 

Corollary 1.33. Let the polynomials p and q be defined in (|1.2[) - (|1.3[l . TTie following are equivalent: 

p 

1) T/ie function R = — is an R-function with exactly n poles, all of which are negative. 

q 

2) The infinite matrix of Hurwitz type H(p,q) defined in (jl.l4[) - (|1.15[) is totally nonnegative and 
Vk+i(p,q) ^ 0. 

3) The finite matrix of Hurwitz type H.k(p,q) defined in (|1.16|) - (|1.17[) is nonsingular and totally non- 
negative. Here k — 2n if degq < degp, and k = In + 1 whenever degg = degp. 

The total nonnegativity of the matrices of Hurwitz matrix also implies the following curious result [16 . 

Theorem 1.34. Let the polynomials p and q defined in (|1.2[) - (11.31) have only nonpositive zeroes, and let 
the function R = qjp be an R-function of negative type. Given any two positive integerij r and I such that 
rl ^ n < (I + \)r, the polynomials 

Pr,i( z ) = o-f)Z l + a r z l ~ x + a 2r z l ~ 2 + . . . + a r i, 
q r ,i(z) = boz' + brz 1 - 1 + b 2r z 1 - 2 + . . . + b rl 

have only negative zeroes, and the function R r j = q r .i/p r .i is an R-function. 

1.4 Stieltjes continued fractions 

Let a real rational function F finite at infinity be expanded into its Laurant series at infinity: 

i^) = s _ 1 + £° + ^ + £a + ... (i.32) 

z z 2 - Z A 



7 The number r is chooser! such that bg + b^. j^ 0. 
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Definition 1.35. The function F is said to have Stieltjes continued fraction expansion if F can be 
represented in the form 



F(z) = c + 



1 



Cj 7^ 0, where T — 



C\Z + 



c 2r if |F(0)| < oo, 

C2r-\Z if -F(O) = oo. 



(1.33) 



C-2 



C 3 Z 



1 

T 



Here r is the number of poles of F, counting multiplicity. 

There exists the following criterion for a rational function to have a Stieltjes continued fraction expan- 
sion [U] (see also [15]). 



Theorem 1.36. Suppose that a rational function F defined in (jl.32[) has exactly r poles. The function F 
has a Stieltjes continued fraction expansion (|1.33p if and only if it satisfies the conditions 



D d (F)^0, i = l,2,...,r, 

S,(F)^0, j = l,2,...,r-l, 

Dj{F) =D j (F) = 0, j = r+l,r + 2,..., 

where Dj (F) and Dj (F) are defined in (| 1 . 5|) and (jl.8p , respectively. 

The coefficients of the continued fraction (|1.33[) can be found by the formules 



(1.34) 

(1.35) 
(1.36) 



c 2j 



C2J-1 = 



D 2 AF) 



Di-i{F) ■ Dj{F) 



, i = i,2 



}***)') 



, j = l,2,...,r, 



(1.37) 
(1.38) 



where D (F) = D (F) = 1. 



We note that the determinant D r {F) in Theorem 1 1 . 361 can be equal zero. But in this case (and only in 
this case), the function F has a pole at according to Corollary II .51 and therefore T = c<z r -\Z in (|1.33l) . 



Let again the polynomials p and q be defined in (II. 2ft (11.31) . Consider the function R = — and suppose 

P ^^ 

that R has a Stieltjes continued fraction expansion (ll.33|) . where r ^ n = degp. Then from (|1.18l) - (|1.21j) 

and (fl~37l) - pT38l) we obtain: 



if degq(z) < degp(z), then 

C-y-i = 
C 2 j = 

if degq(z) = degp(z), then 

C2j-1 = 

C 2 j = 



vi,-i(p,<i) 



V2j-2{p,q) "ri2j(p,q)' 



vh(p,<i) 



r)2j-i(p,q)-r)2j(p,q) 



i = 1,2,. 

J = 1,2,. 



(1.39) 
(1.40) 



vlj(p,q) 

r) 2j -i(p,q) -ri2]+i{P,qy 



»?2i+l (Pi ?) 



ii2j{p,q)-'n2j+2(p,q) ) 



.7 = 1,2,..., 

J=0,l,2,.. 



(1.41) 
(1.42) 
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where k = 2r — 1 if i?(0) = oo, k — 2r if |i?(0)| < oo, and rjj(p, q) are the leading principal minors of the 
infinite Hurwitz matrix H(p,q) (see Definition II .13[) . Here we set rjo(p,q) = 1, and [p] denotes the largest 
integer not exceeding p. 

From Theorems 1 1 . 181 \1 . 241 and 1 1.361 and from the formulae f|1.37[) — f|1.38[) one obtain the following theo- 
rem [IB] , 

Theorem 1.37. Let the function R with exactly r poles, counting multiplicities, be defined (|l.ip . If R has 
a Stieltjes continued fraction expansion (|1.33l) , then R is an R- function if and only if 

c 2j --i>0, j = l,2,...,r. (1.43) 

Moreover, the number of negative poles of the function R equals the number of positive coefficients c 2 j , 
j = 1,2, ... , k, where k = r, if |i?(0)| < oo, and k = r — 1, if R(0) = oo. 

Note that every i?-function with poles of the same sign always has a Stieltjes continued fraction ex- 
pansion. 

Corollary 1.38 (A. Markov, Stieltjes, [27, 38, 39 ). A real rational function R with exactly r poles, counting 
multiplicities, is an R-function with all nonpositive poles if and only if R has a Stieltjes continued fraction 
expansion (|1.33[) . where 

c % > 0, i= l,2,...,2r-l. 

Moreover, if |-R(0)| < oo, then C2- r > 0. 

Corollary 1.39 ([16 ). A real rational function R with exactly r poles, counting multiplicities, is an R- 
function with all nonnegative poles if and only if R has a Stieltjes continued fraction expansion (|1.33l) , 
where 

(-l)'- 1 c i >0, i = l,2,...,2r-l, 

and if |i?(0)| < oo, then C2 r < 0. 

2 Even and odd parts of polynomials. Associated function 

Consider a real polynomial 

p(z) = a z n + a 1 z n ~ 1 H \- a n , Oi, . . . , a n € M, ao > 0. (2.1) 

In the rest of the paper we use the following notation 

1 = 



(2.2) 



where n = degp, and [p] denotes the largest integer not exceeding p. 
The polynomial p can always be represented as follows 



p{z)=p Q {z 2 )+zp 1 (z 2 ), (2.3) 



where 

for n = 2L 



and for n = 21 + 1, 



Po(u) = aou 1 + a 2 u l l + . . . + a„, 
pi(u) = aiu 1 ^ 1 + a3U l ~ 2 + . . . + a„- 

Po(u) = aiu 1 + a^u 1 ^ 1 + . . . + a„, 
Pi(u) = a u l + a 2 u l ~ 1 + . . . + a„_] 



(2.4) 



(2.5) 
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The polynomials Po(z 2 ) and Pi(z 2 ) satisfy the following equalities: 

p(z)+p(-z) 



Po(z 2 ) 

Pi(z 2 ) 



p(z) -p(-z) 



Introduce the following function^: 



$(u) = 



2-- 

Piju) 
Po{u)' 



(2.6) 



(2.7) 



Definition 2.1. We call $ the function associated with the polynomial p. 
From (|2.6[) and ()2.7[) one can derive the following relations: 

P(~z) 

P( z )-p(- z ) 



1 



z<P{z A ) 



P( z )+P(- Z ) 



1 + 



p( z ) 
Pi-z) 

p(z) 



(2.8) 



p{— z) 1 - z<J>(z 2 ) 



p(z) 1 + z<J>(z 2 ) ' 

Now let us introduce two Hurwitz matrices associated with the polynomial p. 
Definition 2.2. The matrix 

/a a 2 a 4 a 6 a 8 ai x 

a\ a^ a§ a-j ag 



Hooip) 



ao a 2 a4 a§ as 
ai a3 a§ 07 



V 



(2.9) 



is called the infinite Hurwitz matrix of the polynomial p. The leading principal minors of the matrix 
H ao (p) will be denoted by T]j(p), j = 1, 2, . . . 

Remark 2.3. According to Definitions 11.131 and 12.21 we have -ffoo(p) = H(po,pi), where po and Pi are the 
even and odd parts of the polynomial p, respectively. 

Together with the infinite matrix H^p), we consider its specific finite submatrix. 
Definition 2.4. The n x n matrix 

(a\ 03 05 CI7 \ 



ft„(p) = 



a 


« 2 


«4 


ttfi 





«1 


03 


a 5 





«() 


a 2 


0.4 














V° 


















fln-2 ««/ 



O-n-1 



(2.10) 



is called the finite Hurwitz matrix or the Hurwitz matrix. 

The leading principal minors of this matrix we denote by A j (p) : 



A» 



where we set a; = for i > n. 



ai a 3 a 5 a 7 

a a 2 a 4 a 6 

a\ a$ a$ 

ao a 2 a 4 







. • a 2 j- 


-1 


. • a 2 j- 


-2 


. • a 2 j- 


-3 


. • a 2i _ 


-4 


a 3 





(2.11) 



3 In the book |12| Chapter XV], F. Gantmacher used the function 



Pl{-u) 

P0 (-«)' 
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Definition 2.5. The determinants A,(p), j = 1, . . . ,n, are called the Hurwitz determinants or the Hurwitz 
minors of the polynomial p. 

Remark 2.6. From Definitions 11.151 and 12.41 it follows that H n (p) = H 2 i(po,pi) if ft = 2/, and H n {p) = 

H 2 i+i(po,Pi) if n = 21 + 1, where po and pi are defined in (|2.3[) — (|2.4|) . 

Obviously, 

Tij(j>) = aoAj-ifr), j = l,...,n, (2.12) 

where Ao(p) = 1. 

Suppose that degpo ^ degpi and expand the function $ into its Laurent series at infinity: 

$ (") = —r-; = s -i + — + — + — + — + •••' ( 2 - 13 ) 

Po(uj u w u u 

where s_i ^ if degpo = degpi, and s_i = if degpo > degpi. 

From (|1.18|l - (|1.23|l we obtain the following relations [HI [12] between the determinants Dj(<&), Dj(<f>) 
defined in (| 1 . 5|) and (jl.8p . the Hurwitz minors Aj(p), and the determinants rjj{p) (J = 1, 2, . . . , n). 

1) If n = 21, then 

A 2j _ 1 (p) = ao 1 V2 j (j>) = % j - 1 D j ($), . 

, - J ^ 1, 2, . . . , (; (■'•14J 

Ajy(p) = e^ Wifr) = (-lYagDjW, 

2) Ifn = 2£ + 1, then 

A 2j (p)=a - 1 r ?2j+1 (p)= f^J £,-(*), j = l,2,...,Z; 
A 23+1 (p)=a^ V23+2 (p) = (-iy^y £,•(*), j = 0,l,...,Z, 



(2.15) 



where A)($) = 1- 
By (|23)) - ([2~5|) and (|2T7)) and by Theorem O and Corollary HH we have 

D j ($) = D j ($) = 0, j>J, (2.16) 

where I is defined in (12.21) . So in the sequel, we deal only with the determinants Dj($), Dj($) of order at 
most I. 

Also in the sequel, we deal only with the determinants r]j(p) of order at most n + 1, since by the 
formulae ([TT5]) - ([OT]l and ([2~TB|l and by Remark [231 we have 

n ] (p)=0, for j>n+l. (2.17) 

In Section [S]we also consider the case when n = 21 + 1 with a! = and 03 7^ 0. In this case, we have 
degpo = dcgp! — 1 = I — 1, so the function $ has the form 

*/ \ Pi( u ) , , so , Si s 2 s 3 . no"' + azu^ 1 + . . . + a 2l 

®( u ) = — tt = s -2 u + s-i + — + — + — + — + ••■ = — rrr — r=rz — z • ( 2 - 18 ) 

Po( u ) u u u A u 4 0,311' i + a^u 1 z + . . . + 0,21+1 



_ . , ao , _, . ao P2(w) , 

It is easy to see that t- 2 = — , so we can represent $ as a sum <P{u) = — u-\ — — -, where 

a 3 a 3 po(u) 

1 \ ( «o \ /—1 / ao A 1-2 ( a Q \ 

p 2 (u) = [a 2 a 5 ]u + a 4 a 7 )u +...+ [ a 2i - 2 a 2 /+i u + a 2 i- 

\ a 3 J \ a 3 J \ a 3 J 

Note that for the function 

rv \ P2 ^ if,/ ^ a ° , s ° 1 Sl 1 S2 , S3 , 

8« := — rr = $ M ) u = S-i-\ h -j + -~ + -7 + • • • , 

PoW 03 u w u m 4 
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the following equalities hold 

D j (Q)=D j ($), and Dj (6) = Dj ($) for j = 1,2, 



(2.19) 



By Definition 11.151 and by the formulas (|1.20[) . and (jl.231) applied to the polynomials p and p 2) we have, 
for j = 1,2,..., 



,2j+l ; 



-a a 3 V2j 


bo 


P2) = - 


-a a 3 


02 


03 
a 

05 

03 






a 5 
a 

04 07 

03 
a 3 
a 

a 2 a 5 ... 

03 




04j + l 

a 

04j O4J+3 

«3 

04J-1 

ao 

04J-2 a4j+i 

03 

a 

a 2j+2 02J+5 

a 3 


a 3 




« 5 


a 7 




a 4j+3 




ao 02 




0,4 


«(i 




04J+2 










a.2 


«3 

a 

a 5 

03 





a 5 
a 

04 07 

03 

a 3 
a 
a 2 a 5 ... 

03 


04j + l 

a 

04j 04j+3 

«3 
O4J-I 

a 
04j-2 a4j + i 

03 


= 

















a 

02J+2 0-2J+5 

«3 





u 


03 


05 


07 


«o 


02 


04 


o 6 








03 


«5 





«o 


«2 


0,4 











«3 








«o 


«2 















(p) 



04.7+3 

04 j +2 
04J+1 

a 4j - 

a 4j _l = A 2;-2 

a 4j-2 
02j+2 

where we set a, = for i > 21 + 1. 

Denote the coefficients of the polynomial p 2 by bj : 

Oj = a 2j+ 2 a 2 j+5, J = U, 1, . . . 

«3 

where 022+3 — 0. Then the function u$(u) has the form 

Kf \ "o 2 , & o . Pa( u ) 

a 3 a 3 po« 



'OoOg' 



2J ' +1 ^($), 



,1-1, 



(2.20) 



b .._ ,.. s /'r fy) _ \ . ;_i , A_ ^ 



where 

p 3 (u) := up 2 (u) 
By Definition 11.151 and by the formulas (|2.19p , (|1.20p and (|1.23p applied to the polynomials po an d _p 3 we 



up (u) = [bi a 5 m + fr 2 a 7 u + 

03 \ a 3 / V a 3 



fr(-l 02i+l u. 

«3 



get (putting bi = for i > I — 1) 
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aoa 3 



2j+2 D j (Q) = a a 3 ! V(po,P3) 



aoa 2 . 



a a 3 



«3 









as 

«3 

6o 



«7 
«5 

bo 



L u L u 

Ol 05 02 07 



«3 







a3 
a3 

oi a 5 

a 3 



= (-1) J+1 



U a 3 
a a 2 




04 



«3 

Oi as 

a 3 










dij+3 

a 4j+l 

t>2j a 4j+3 

0-3 

aij-i 

h b ° 

02j-i - — a 4 i+i 

«3 



h b ° 

O7+I a 2j+5 

«3 

a 7 

a 6 



62 


a 5 

60 

a 7 .. 

a 3 


h b ° 

t>2j a 4j"+3 

a 3 


01 


a 3 

bo 

a 5 . . 

03 


04.7-1 

02J-1 a 4j+l 

«3 



b 7+ i - 



«3 



■ 02.7+5 



= (-l) J a a 3 



«3 


fl 5 


«7 


00 


01 


&2 





«3 


a 5 





00 


01 








«3 








00 









02 


( 















































«3 
«0 



a-, 
a 

":j 03 

03 
ao 
«2 «5 04 

03 



a 2 



«7 

a 
09 

03 
a 5 

a 
a 7 

03 

03 
a 

o 5 

03 









04.7+5 
04J+4 

04.J+3 
Oo 

04j+2 04j+5 

03 

O4.7+I 
O 

04j 04^+3 

03 
04j— 1 

a 

O4J-2 04j+i 

03 



02j+3 



04.7+3 

b 2 j 
O4J+I 
623-I 
04j_l 

&2.7-2 



°2j+3 



(-iy+ 1 A 2j +3(p) = aoa^ J D J (ci>), 



where we set cjj = for i > 2/ + 1. Here we used the equalities Dj(Q) = Dj(uQ), j = 1, 2, . 
Thus, we obtain the following formulae 



AJ«(P) = o ^+iO) = -o a5 J ^j-iC^) = «3 J z A 2 (p)^--i($), 

A 2j+ i(p) = a^ V2j+2 (p) = (-lyaoafD^) = (-l^o^AafcO^-iC*), 

where r?i(p) are the leading principal minors of the matrix H oc (p) defined in 
formulae 

A 2 (p) = Oq 1 r} 3 (p) = -a a 3 , A 3 (p) = a 1 7y 4 (p) = -a a a\ 

that follow from the equality a\ = 0. 

Note that, as above, the following holds: 

DjW = £,-(#) = 0, j>/-1. 



.7 = 2,3,..., /, 



(2.21) 



.7 = 2,3,...,?, 

Here we used the 

(2.22) 



3 Stable polynomials 

This section is devoted to some basic facts of the theory of Hurwitz stable and quasi-stable polynomi- 
als. We expound those facts from the viewpoint of the theory of i?-functions. 
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3.1 Hurwitz polynomials 

Definition 3.1 ([T71IH]). The polynomial p defined in (j2.ll) is called Hurwitz or Hurwitz stable if all its 
zeroes lie in the open left half-plane of the complex plane. 

At first, we recall the following simple necessary condition for polynomials to be Hurwitz stable. This 
condition is called usually Stodola theorem |12[ [B] . 

Theorem 3.2 (Stodola). If the polynomial p is Hurwitz stable, then all its coefficients are positivqj. 

Now we show a connection between Hurwitz polynomials and i?-functions. This connection will allow us 
to apply all statements from the Section Q] to Hurwitz polynomials to obtain such basic criteria of Hurwitz 
stability as Hurwitz criterion and Lienard and Chipart criterion (Theorems (I3.5[) and 13.71 respectively). 
We prove the following theorem using a method suggested by Yu.S. Barkovsky [5 . 

Theorem 3.3 f |121 15]). A real polynomial p defined in (|2.ip is Hurwitz stable if and only if its associated 
function $ defined in (|2.7[) is an R-function with exactly I poles, all of which are negative, and the limit 



lim $(u) is positive whenever n = 21 + 1. The number I is defined in (|2 .2|) 

u— »±oo 

Proof. Let the polynomial p be Hurwitz stable, that is, 

p(A) = =>> RcA<0. 
First, we show that 

< 1, \/z : Rez > 0. 



p(-z) 



p{z) 



(3.1) 



(3.2) 



Note that the polynomials p(z) and p(— z) have no common zeroes if p is Hurwitz stable, so the function 

— — — has exactly n poles. The Hurwitz stable polynomial p can be represented in the form 

p(z) 



p(z) =a Y[(z- A fc ) Y[(' 



•&)(* 



where A& < 0, Re£j < and Im£j ^ 0. Then we have 

p(-z) 



p(z) 



It is easy to see that the function of type 
plane to the unit disk. In fact, 

2 



k ' 

z + a 
z — a 



Afc|-|-rl* + 6lF + & 



Afei -*•:*■ \z 



(3.3) 



•£ill*-&l 
where Re a < 0, maps the right half-plane of the complex 



z + a 



a 



(Re z + Rca) + (Im z + Im a) 
(Re z — Re a) 2 + (Im z + Im a) 



< 1 



whenever Re 2: > and Re a < 0. 



u( — z] 
Now from (I3.3[) it follows that the function — -—^ also maps the right half-plane to the unit disk as 

P[ z ) 
a product of functions of such a type. Thus, the inequality Q3.2p is valid. 

1 — z 

conformally maps the unit disk to 



At the same time, the fractional linear transformation z n- 
the right half-plane: 

\z\ < 1 =$> Re 



1-z 



1 + 2 
1-UI 2 



1 + Zj |1 + 2| 2 >0 - ^ 

Consequently, from the relations (|2.8|) . p.2| and p.4| we obtain that the function z$(z 2 ) maps the right 
half-plane to itself, so the function — z$(— z 2 ) maps the upper half-plane of the complex plane to the lower 
half-plane: 



Imz > 



Re(-iz) > 



Re [-iz<5>{-z 2 )} = Im [z${-z 2 )} > =*• Im [-z$(-z 2 )] < 0. 



3 More exactly, the coefficients must be of the same sign, but ao > by 12.lt . 
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Since p is Hurwitz stable by assumption, the polynomials p(z) and p(— z) have no common zeroes, therefore, 
Po and pi also have no common zeroes, and Po(0) ^ by (|2.3p . Moreover, by Theorem 13.21 we have a > 
and ai > 0, so degp — I (see (|2.4p - (l2.5p ). Thus, the number of poles of the function — z&(— z 2 ) equals 
the number of zeroes of the polynomial po(— z 2 ), i.e. exactly 21. 

So according to Theorem 1 1.1 81 the function — z$(— z 2 ) can be represented in the form (jl.251) . where all 
poles are located symmetrically with respect to and /3 = 0, since — z$(— z 2 ) is an odd function. Denote 
the poles of — z$(— z 2 ) by ±v\, . . . , ±i/i such that 

< Vy < V<z < . . . < Vi . 

Note that v\ ^ 0, since po(0) / as we mentioned above. 

Thus, the function — z$(— z 2 ) can be represented in the following form 

ii i „ 

-z$(-z 2 ) = -az + V— ^ hV 7j = -az + V 9 7jZ ? , a ^ 0, jj,Vj > 0. 

v 7 ^ z - v„ ^ z + Vl ^ z 2 - v 2 u J 

Dividing this equality by —z and changing variables as follows — z 2 — > u, 27j — > /3j, in — > wj, we obtain 
the following representation of the function $: 



where a ^ 0, /3, > and 



$ (u) = ^ = a + V^^, (3.5) 

po(«) £j« + <"i 



< Wi < U>2 < ■ ■ ■ < Wj . 

Here a = whenever ?i = 2^, and a = — > whenever n = 21 + 1. Since $ can be represented in the 

ai 

form (13. 5|) . we have that by Theorem 1 1.1 8[ $ is an i?-function with exactly I poles, all which are negative, 

and lim $(u) = a > as n = 21 + 1. 

u— >±oo 

Conversely, let the polynomial p be defined in (|2.1j) and let its associated function $ be an i?-function 
with exactly I poles, all of which are negative, and lim <fr(w) > as n = 21 + 1. We will show that p is 

u— >±oo 

Hurwitz stable. 

By Thcorcm ll.181 $ can be represented in the form (13.51) . where a = lim <&(u) Js such that a > 

if n = 2Z + 1, and a — if n — 21. Thus, the polynomial po has only negative zeroes, and the polynomials 
Po and P\ have no common zeroes. Together with (|2.3p and (|2.7p . this implies that the set of zeroes of the 
polynomial p coincides with the set of roots of the equation 

z$(z 2 ) = -l. (3.6) 

Let A be a zero of the polynomial p and therefore, of the equation (|3.6p . Then from (|3.5p and (|3.6p we 
obtain 

|2 i , , 



-1 = Re [A$(A 2 )] = 



■£* 



. L J |A2 + Wi |a 



Re A, 



where a ^ 0, and /3j, ujj > for j = 1, . . . , /. Thus, if A is a zero of p, then Re A < 0, so p is Hurwitz 
stable. □ 

Consider the Hankel matrix S(§) = ||si+j||^- =0 constructed with the coefficients of the series (|2. 13|) 

and its determinants Dj($>) and Dj($>) defined in (jl.51) and (ll.8p . respectively. From Theorems II . 181 [3~3l 
and 11.121 and from Corollary 11.261 we obtain the following Hurwitz stability criteria: 



Theorem 3.4. Let a real polynomial p be defined by (|2.ip . The following conditions are equivalent: 
1) the polynomial p is Hurwitz stable; 
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2) the following hold 



s-i > for n = 2l + l, 




D j (<b)>0, j = l,...,l, 




(-l)'Z>,-(fc) >0, j = l,.. 


.,1, 



(3.7) 



where I 



3) the matrix S(Q) is sign regular of rank I, and s_i > for n = 21 + 1. 

Proof. In fact, by Theorem l3.3l the polynomial p is Hurwitz stable if and only if its associated function $ 
is an i?-function with exactly I poles, all of which are negative and lim $(u) = s_i ^ 0. According to 

u— >±oo 

Theorem 11.181 and Corollary 11.261 this is equivalent to the inequalities (|3.T[) . But these inequalities also 
are equivalent to the sign regularity of the matrix 5($) by Theorem II. 121 □ 

Our next theorem provides stability criteria in terms of coefficients of the polynomial p. 

Theorem 3.5. Given a real polynomial p of degree n as in (j2.ll) . the following conditions are equivalent: 

1) the polynomial p is Hurwitz stable: 

2) all Hurwitz determinants Aj(p) are positive: 

Ai(p) > 0, A a (p) > 0, . . ., A n (p) > 0; (3.8) 

3) the determinants rfj (j>) are positive up to order n + 1 : 

m { P ) > 0, m {p) > 0, . . . , Vn+1 {p) > 0; (3.9) 

4) the matrix H n {p) defined in (|2.10[) is nonsingular and totally nonnegative; 

5) the matrix H^p) defined in (|2.9[) is totally nonnegative with nonzero minor r) n+ i(p). 

Note that the equivalence of 1) and 2) is the famous Hurwitz criterion of stability. The implications 
1) => 4) and 1) => 5) were proved in [3J[3D]. The implication 4) =>• 1) was, in fact, proved in [3J. 
However, the implication 5) => 1) is probably new. 



Proof. Indeed, by Theorem 13.41 the polynomial p is Hurwitz stable if and only if the inequalities (|3.7p 
hold. According to the formulas (|2.14[) - ([2.15[) . these inequalities are equivalent to (|3.8p . By (|2.12[) . the 
inequalities (|3.8p are equivalent to (|3.9p . since r]i(p) = ao > (sec (|1.14l) - (|1.15p and Remark l2.3p . 

Furthermore, by Theorem l3.31 the polynomial p is Hurwitz stable if and only if its associated function $ 
is an i?-function with exactly I = 



21 + 1. According to Corollary 



all of which are negative, and lim <&(u) = — = s_i > if 

u— >±oo a\ 

1.331 and Remark 12.31 this is equivalent to the total nonnegativity of 



the matrix H aD (p). Moreover, r) n +i(p) 7^ by Corollary 11.331 Thus, the condition 1) is equivalent to the 
condition 5). 

The conditions 1) and 4) are equivalent also by Corollary 11.331 and Remark 12.61 □ 

Remark 3.6. Note that total nonnegativity of the Hurwitz matrix of a Hurwitz polynomial was established 
in [3J [201 H • 

Now we are in a position to prove a few another famous criteria of Hurwitz stability, which are known as 
the Lienard and Chipart criterion and its modifications (sec [25, 12J). These criteria are simple consequences 
of Theorems [QD1 IOT1 and [ 



Theorem 3.7. The polynomial p given by (|2.1[) is Hurwitz stable if and only if one of the following 
conditions holds 

1) a n >0,a„_ 2 > 0,a„_ 4 > 0, . . . , A„_i(p) > 0,A n _ 3 (p) > 0,A„_ 5 (p) > 0, . . .; 



20 

2) a n >0,a„_! >0,a„_ 3 > 0, . . . , A n _ x {p) > 0,A„_ 3 (p) >0,A„_ 5 (p) >0,...; 

3) a n >0,a„_ 2 >0,a„_ 4 >0,..., A„(p) > 0, A„_ 2 (p) > 0, A n _ 4 (p) > 0, . . .; 

4) a„ >0,a n _i > 0,a„_ 3 > 0, . . . , A n (p) > 0, A„_ 2 (p) > 0, A„_ 4 (p) > 0, . . . 

Proof. Let n = 21. By Definitions 11.151 and 12.41 and by Remark 12.61 we have H n (p) = H2i(po,pi), where 
the polynomials po an d P\ are defined in (|2.3p - (|2.4j) . Obviously, Aj(p) = Aj(po,pi), j = l,...,n. The 
assertion of the theorem follows now from Theorems 13.31 and 11.301 

In the case n = 21 + 1, the assertion of the theorem can be proved in the same way using Theorem ll.311 
instead of Theorem 11.301 □ 

We now recall a connection between Hurwitz polynomials and Stieltjes continued fractions (see |12j). 

Theorem 3.8. The polynomial p of degree n ^ 1 defined in (|2.ip is Hurwitz stable if and only if its 
associated function $ has the following Stieltjes continued fraction expansion: 

$>(u)=c ^ , with c t > 0, i = l,...,2l, (3.10) 



C\U 



C>2 



C 3 U 



1 



1 

C2l-lU-\ 

C21 



where Co = if n is even, cq > if n is odd, and I as in 

Proof. In fact, by Theorem 13.31 the polynomial p is Hurwitz stable if and only if the function $ can be 
represented in the form (|3.5[) , where a — c$ = if n — 21, and a = cq = — > if n = 21 + 1. Now 
the assertion of the theorem follows from Theorem 11.181 and Corollarv ll.381 □ 

From the formulae (ll.37[) - (|1.38p it follows that the coefficients Ci of (13.10)) can be found by the following 

formulas 

£>?($) 
^=- n '. g— , j = l,2,...,l, (3.11) 

C ^ = D-M^wr i=1 ' 2) •-'' (3 ' 12) 

where D ($) = Z> ($) = 1. If n is odd, then c = — > 0. 

Qi 
Using the formula? (|3.11[) — (13.121) and (|2.14[) — (|2.15l) . we can represent the coefficients Cj in terms of 

Hurwitz determinants Ai(p): 
1) If n = 21, then 

A ^ l(:P) i = l,2,...,n, (3.13) 



2) If n = 21 + 1, then 



A^ 2 (p) ■ AiipY 

A?(p) 
A i _ 1 (p)^A i+1 (J^) , 



i = 0,l,2,...,n, (3.14) 



where we set A_i(p) = — and Ao(p) = 1. 

At last, from Theorems 11.221 and 13.31 one can obtain the following simple result. 
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Theorem 3.9. Let p be a Hurwitz polynomial of degree n ^ 2 as in (|2.1|) . Then all the polynomials 



n-2j 
i=0 



n — i 


( 


n — i 


2 


2 



1. 



also are Hurwitz stable. 

Proof. By Theorem 13.31 if p is Hurwitz stable, then the function $ = pi/po is an R- function. According 
to Theorem 1 1.22[ all functions $j = p\ 3 /pq , j = 1, . . . , — — 1, are i?-functions with negative zeroes and 



n 
2 

D 



poles. Now Theorem 13.31 implies that all polynomials 

Pj (z)=p^(z 2 )+zp[ j \z 2 ), J = l, 
are Hurwitz stable, as required. 



3.2 Quasi-stable polynomials 

In this section we deal with polynomials whose zeroes lie in the closed left half-plane of the com- 
plex plane. 

Definition 3.10. A polynomial p of degree n is called quasi-stable with degeneracy index m, ^J m ^ n, 
if all its zeroes lie in the closed left half-plane of the complex plane and the number of zeroes of p, counting 
multiplicities, on the imaginary axis equals m. We call the number n — m the stability index of the 
polynomial p. 



Throughout this section we use the following notation 



(3.15) 



where n and m are degree and degeneracy index of the polynomial p, respectively. 

Obviously, any Hurwitz polynomial is quasi-stable with zero degeneracy index, that is, it has the 
smallest degeneracy index and the largest stability index (which equals the degree of the polynomial). 
Note that if the degeneracy index m is even, then p(0) ^ 0, and if m is odd, then p must have a zero at 0. 

Moreover, if p is a quasi-stable polynomial, then 



p(z) = p (z 2 ) + z Pl (z 2 ) = f(z 2 )q(z) = f(z 2 ) [q a (z 2 ) + z qi (z 2 )] 
m~\ 
~2 



where f(u) is a real polynomial of degree — with nonpositive zeroes, and q is a Hurwitz stable poly- 
nomial if m is even, and it is a quasi-stable polynomial with degeneracy index 1 if m is odd. Using this 
representation of quasi-stable polynomials, one can extend almost all results of Section l3~T1 to quasi-stable 
polynomials in the same way, so we state them here without proofs. We only should take into account 
that the function 

Pi( u ) f{u)qi{u) gi(u) 



$(m) 



(3.16) 



p (u) f(u)q {u) q (u) 
has a pole at whenever p(0) = 0. 

The analogue of Stodola necessary condition for quasi-stable polynomials is the following. 
Theorem 3.11. If the polynomial p defined in (12.11) is quasi-stable, then all its coefficients are nonnegative. 
The next theorem is an extended version of Theorem [ 



Theorem 3.12. The polynomial p of degree n defined in (j2.ll) is quasi-stable with degeneracy index m if 
and only if the function $ defined in (|3.16p is an R-function of negative type with exactly r poles all of 
which are nonpositive, and lim $(u) is positive if n = 11 + 1. The number r is defined in (|3.15p . 

u— >±oo 



If we expand the function $ into its Laurent series (|2.13p . then the Hankel matrix S{<&) = ||sj+j||o° has 



rank r, where r as in (I3.15[) . 
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Theorem 3.13. Let the polynomial p be defined in (|2.1[) . The following conditions are equivalent: 

1) the polynomial p is quasi-stable with even degeneracy index m; 

2) the following hold 

s_i > for n = 21 + 1, 

£>,($) >0, i = l,...,r, 

(-1^(4) > 0, i = l,...,r-l, 

and (— l) r D r ($) > /or even m, but D r (&) = /or odd m. Here r is defined in (| 3 . 1 5 1) .' 

5j the matrix S(<l>) is sign regular of rank r for even m, S(<&) is r — 1-sign regular of rank r for odd m, 
and S-i is positive if n = 11 + 1. 



It is also easy to extend Theorem 13.51 to quasi-stable polynomials. 

Theorem 3.14. Let the polynomial p of degree n be defined in (|2.1[) . The following conditions are equiv- 
alent: 

1) the polynomial p is quasi-stable with degeneracy index m; 

2) determinants Aj(p) are positive up to order n — m: 

Ai(p) > 0, A 2 (p) > 0, . . . , A„_ m (p) > 0, A„_ m+1 (p) = . . . = A„(p) = 0; 

3) determinants rjj (p) are positive up to order n — m + 1 : 

rji{p) > 0, n 2 {p) > 0, . . . , 77„_ m+ i(j>) > 0, r? n _ m+i (p) =0, i = 2, 3, . . . ; 

^ i/ie matrix if oo (p) is totally nonnegative and 

Vn- m +l{p) ¥" 0, Vn-m+i(p) = 0, i = 2, 3, . . . 

The implication 1) =^> 4) was proved, indeed, in [31 [50] • The implication 4) =>■ 1) seems to be new. 
Thus, we established that a polynomial is quasi-stable if and only if its infinite Hurwitz matrix is totally 
nonnegative. 

The following interesting property of quasi-stable polynomials is a simple consequence of Theorems ll.341 
and 13.31 



Theorem 3.15. Let the polynomials p of degree n — 21 defined in (|2.1[) be quasi-stable. Given any positive 
integer r(^ n), the polynomial 



(\ Ak Ak I A k' 2. 2.K "> A k 4 

z) = a Q z + a 2r -iz + a 2r z + a± r -\z +a ir z + ■ ■ ■ + a 2r k-iz + a 2r k, 



where k = 



, also is quasi-stable. 



Theorem 3.16. Let the polynomials p of degree n = 21 + 1 defined in (|2.1[) be quasi-stable. Given any 
positive integer r(^ n), the polynomial 

p r (z) = a Q z 2k+1 + aiz 2k + a 2r z 2k ~ 1 + a 2r +iz 2k ~ 2 + a± r z 2k ~ 3 + a ir+ iz 2k ^ 4 + . . . + a 2rk z + a 2rk+ i, 



where k 



, also is quasi-stable. 
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Although total nonnegativity of the infinite Hurwitz matrix is equivalent to quasi- stability of polyno- 
mials, total nonnegativity of the finite singular Hurwitz matrix is not equivalent to quasi-stability as was 
noticed by Asner [3 . It is clear from Theorem 13.141 that the finite Hurwitz matrix of a quasi-stable polyno- 
mial is totally nonnegative. However, given a real polynomial p of degree n, if TL n (p) is totally nonnegative, 
p is not undertaken to be quasi-stable. In fact |16) . p can be represented in the form p(z) = q(z)g(z 2 ), 
where the polynomial q is Hurwitz stable, and the polynomial g is chosen such that the matrix Tt n (p) 
is totally nonnegative but g(z 2 ) has zeroes in the open right half-plane. Such choice is possible as was 
mentioned in [16] . 

The connection between quasi-stable polynomials and Stieltjes continued fractions is similar to Hurwitz 
stable polynomials' one. Namely, one can prove the following extended version of Theorem] 



Theorem 3.17. The polynomial p of degree n defined in (|2.1I) is quasi-stable with degeneracy index m if 
and only if the function $> has a Stieltjes continued fraction expansion: 

• («) = c„ + , With c t >0, T={ °»> tf m %S eVen > (3.17) 



C\U 



C 2r _ lU, if m is odd. 



C-2 



C 3 U 



1 

T 



Here Co = if n is even, cq > if n is odd, and r is defined in (|3.15|1 . 

The coefficients a in (|3.17|) also can be found by the formulae (|3.13[l - (|3.14|l . 

Finally, we note that a quasi-stable polynomial p with degeneracy index 1 can be represented as 
follows: p(z) — zq(z), where q is a Hurwitz stable polynomial. Therefore, for quasi-stable polynomials 
with degeneracy index 1, one can establish an analogue of Lienard and Chipart criterion. 

Theorem 3.18. The polynomial p given by (|2.1[) is quasi-stable with degeneracy index 1 if and only if one 
of the following conditions holds 

1) a n =0,a n _i >0,a„_ 2 > 0, a„_ 4 > 0, . . . , Ai(p) > 0,A 3 (p) > 0,A 5 (p) > 0, . . .; 

2) a„ =0,a n _i >0,a„_ 3 > 0,a„_ 5 > 0, . . . , A x (p) > 0,A 3 (p) > 0,A 5 (p) > 0, . . .; 

3) a„ =0,a n _i >0,a„_ 2 > 0, a„_ 4 > 0, . . . , A 2 (p) > 0,A 4 (p) > 0,A 6 (p) > 0, . . .; 

4) a„ = 0, a„„x > 0, a„_ 3 > 0, a„_ 5 > 0, . . . , A 2 (p) > 0, A 4 (p) > 0, A 6 (p) > 0, . . . 

4 Self-interlacing polynomials 

In Section [3] we established that the function $ (see (|2.7[l ) associated with a Hurwitz stable polynomial 
(quasi-stable polynomial) maps the upper half-plane of the complex plane to the lower half-plane and pos- 
sesses only negative (nonpositive) poles. Now we are in a position to describe polynomials whose associated 
function $ also maps the upper half-plane to the lower half-plane but has only positive (nonnegative) poles. 

4.1 General theory 

Definition 4.1. A real polynomial p(z) is called self-interlacing if all its zeroes are real and simple and 
interlace zeroes of the polynomial p{—z). 

In other words, if A, are the zeroes of a self-interlacing polynomial p, then one of the following holds: 

< Ai < -A 2 < A 3 < . . . < (-1)"- 1 A„, (4.1) 

0<-Ai<A 2 <-A 3 <...<(-l) n A„, (4.2) 

where n — degp. 
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Definition 4.2. A real polynomial p of degree n is called self-interlacing of type I (of type II) if all its 
zeroes are real and simple and satisfy the inequalities (|4.1[) (respectively, (J4.2J) 1 ). 

If a polynomial p of degree n is self-interlacing of type I, then its minimal absolute value zero is positive. 
Moreover, if n = 21, then p has exactly I negative zeroes and exactly I positive zeroes. Its zero A„, which 
has the maximal absolute value is negative. If n = 11 + 1, then the polynomial p has exactly I negative 
zeroes and I + 1 positive zeroes. In this case, A„ is positive. 

Note that a polynomial p(z) is self-interlacing of type I if and only if the polynomial p(— z) is self- 
interlacing of type II. So in the sequel, we deal only with self-interlacing polynomials of type I. 

Theorem 4.3. Let p be a real polynomial of degree n ^ 1 as in (|2.ip . The polynomial p is self-interlacing 
of type I if and only if its associated function $ defined in (J2.7P is an R-function with exactly I poles, all 
of which are positive, and lim $(u) is negative whenever n = 21 + 1. 

u— >±oo 

Proof. Let p be self-interlacing of type I. First, we show that the function 

p(-z) 



G(z) = 



p(z) 



is an i?-function. In fact, by Definition 14.11 the zeroes of the polynomials p(z) and p(—z) are real, 
simple and interlacing, that is, between any two consecutive zeroes of one polynomial there lies exactly 
one zero, counting multiplicity, of the other polynomial, so G(z) or —G(z) is an i?-function according 
to Theorem 11.181 By Corollary 11.201 G is monotone between its poles. So it remains to prove that the 
function G is decreasing between its poles. 

Let n = 21, then by (I4.1[) . the maximal pole of G is A„_i > 0, but its maximal zero is — A„ > 0, 
which is greater than A„_i according to Definition 14.21 fsee (I4.1[l ). Therefore, in the interval (— A„,+oo), 
the function G has no poles and zeroes. At the same time, lim G(z) — —1. Consequently, in the 

z— >±oo 

interval (— A n , +oo), G(z) decreases from to —1. So, G is an i?-function. In the same way, one can prove 
that if n = 21 + 1, then G is also an R- function. 

Thus, if p is self-interlacing of type I, then G maps the upper half-plane to the lower half-plane. 
From (|2.8p we obtain that 

-ev£§$ («) 

1 + z 
Since the fractional linear transformation z t— > conformally maps the lower half-plane to the lower 

1 — z 
half-plane: 

fl-z\ 2Imz . , 

Im,<0^Im(^ TT -j^^ T - F <0, (4.4) 

from (I4.3I) - (|4.4I) it follows that the function z$(z 2 ) maps the upper half-plane to the lower half-plane, that 
is, z$(z 2 ) is an i?-function. 

Since p is self-interlacing of type I by assumption, the polynomials p(z) and p(—z) have no common 
zeroes, therefore po and p\ also have no common zeroes and po(0) ^ by (|2.3[) . Therefore, the number of 
poles of the function z$(z 2 ) equals the number of zeroes of the polynomial po(z 2 ). If n — 21, then by (|2.4p , 
degpo = I, so z$(z 2 ) has exactly 21 poles. If n = 21 + 1, then by (J2.5I) . degpi = / and degpo ^ I, so z$(z 2 ) 
has at most 21 poles, and it has exactly 21 + 1 zeroes, since po(0) ^ as we mentioned above. But z$(z 2 ) 
is an i?-function, therefore, it has exactly 21 poles by Theorem 1 1.1 81 

Thus, Theorem 11.181 implies that the function z$(z 2 ) can be represented in the form (|1.25p . where all 
poles are located symmetrically with respect to and (3 = 0, since z$(z 2 ) is an odd function. Denote the 
poles of z$(z 2 ) by ±i/i, . . . , ±i/; such that 

< v\ < v<z < . . . < vi . 

Note that fi^O since po(0) ^ 0. 

So the function z$(z 2 ) can be represented in the following form 

l i ' 2 z 

z$(z 2 ) = -az + V -^— + V -^— = -az + V „ ljZ „ , a > 0, 7,-, Vi > 0. 

j=l J j=l ^ ,=1 3 
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Divide this equality by z and changing variables as follows z — \ u, 2jj —> f3j, v- —* loj, we obtain the 
following representation of the function $: 



l 

L 



$( u ) = -a + V^^, (4.5) 



where a ^ 0, /3j > and 

< Wi < u>2 < . . • < UJi . 

do 
Here a = whenever n = 21, and — a = — < whenever n = 21 + 1. Since $ can be represented in the 

a\ 
form (|4.5p . by Theorem II. 181 $ is an i?-function with exactly I poles. Moreover, from (|4.5[) it also follows 

that all poles $ are positive and lim <J>(u) < if n = 21 + 1. 

u— >±oo 



Conversely, let the polynomial p be defined in (12.1[) . and let its associated function $ be an R- function 
with exactly I poles, all of which are positive, and let lim <f>(w) < when n = 21 + 1. We will show that 

u— >-±oo 

p is self-interlacing of type I. 

By Theorem 1 1.181 $ can be represented in the form (|4.5p . where —a — lim $(w) ^ such that a > 

u— >±oo 

if n = 2Z + 1, and a = if n = 21. Thus, the polynomial p nas on ly positive zeroes, and the polynomials 
Po and V\ have no common zeroes. Together with (I2.3P and (J2.7I) . this means that the set of zeroes of 
the polynomial p coincides with the set of roots of the equation 

z$(z 2 ) = -1. (4.6) 

At first, we study real zeroes of this equation, so we consider only real z. Since pq has only positive zeroes, 
we have p(0) = po(0) ^ 0. Thus, we put z € K\{0}. Changing variables as follows z 2 — > u > 0, we rewrite 
the equation (|4.6|) in the following form, which is equivalent to (|4.6[) for real nonzero z: 




u>0. (4.7) 



Note that all positive roots (if any) of the first equation (14.71) are squares of the positive roots of the equa- 
tion (|4.6[) , and all positive roots (if any) of the second equation (|4.7p are squares of the negative roots of 
the equation (|4.6p . 

Let n — 21. Then lim $(u) = —a = 0. Consider the function Fi(u) = $(u) -\ -= whose 

set of zeroes coincides with the set of roots of the first equation (14.71) . The function F\{u) has the 
same positive poles as $(w) does, that is, wi,W2, ■ • ■ , wj. Moreover, F\{u) is decreasing on the inter- 
vals (0,wi), (uj\,u)2), ■ • ■ j (wj— ii wz), (i^ 1 ;, +oo) as a sum of functions that are decreasing on those intervals. 
Since F\ is decreasing on (wi,+oo) and F\(u) — > +0 as u — > +oo, we have -Fi(w) > for u € (oj;,+oo). 
Further, Pi(u) — >• — oo asu/ w; and F\{u) — > +oo as u \ Wj, i = 1, . . . , I. Also F\{u) — > +oo when- 

< oo and —= — > +oo as u — > +0. Thus, on each of the intervals 



Pi(0) 



ever u — > +0, since |$(0)| 

Po(0) 

(0,Wi), (oj 1 ,aj 2 ), ■ • • , (w;_i,o;;), the function i<\ decreases from +oo to — oo, and therefore, it has exactly 
one zero, counting multiplicity, on each of those intervals. Denote the zero of F± on the interval (wj_i, Wj) 
by /if, i — 2, . . . ,1 (/i.j > 0). Also we denote by /if (/ii > 0) the zero of F\ on the interval (0, Wi). Since 
$ is an .R-function by assumption, it has exactly one zero, counting multiplicity, say &, on each of the 

intervals (wj, Wj+i), i = 1, . . . ,1 — 1. But i^ife) = — - f= > 0, so we have 

n\ < wi < Ci < /"I < w 2 < & < • • • < 6-i < M? < w«- ( 4 - 8 ) 

Thus, the first equation (J4.7I) has exactly Z positive roots fa, all of which are simple and satisfy the 
inequalities 
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Consider now the function -FV - ") = $(u) p= whose set of zeroes coincides with the set of roots of 

1 Pi 
the second equation (14. 711 . Since $(0) = — > — < by (I4.5I) . and $(u) ->-oo as«/wi, the function 

z— 1 

$ decreases from $(0) < to — oo on the interval (0,wi). Consequently, ^(u) < for all u € (0,wi). It 
is clear that F2OU) — > +00 as « \ w, and ^(u) — )■ —00 as u /* u>j, ? = 1, . . . , ?. Since i^f/uf ) = < 0, 

Mi 

the function F2 has an odd number of zeroes, counting multiplicities, on each of the intervals (wj,//|), 

2 

i = 1, . . . ,1 — 1. Besides, F 2 (m) = Fi(u) = < whenever u e \^f,oj i+ i), i = 1, . . . ,1 — 1, because 

V u 

Fi(u) < on those intervals. More exactly, since -F^C^i) = —j= < 0, * = 1, . . . , r— 1, the function -F^u) has 

an odd number of zeroes on each interval (uii, £,i) and an even number of zeroes on each interval [£j, Aif +1 ), 
i = 1, . . . , r — 1. So F2 has at least I — 1 zeroes in the interval (wi, Wj). Consider now the interval (wj, +00) 
and show that i^Ci-O — > — • as u — > +00. In fact, since the function $ is an i?-function by assumption, 
we have Sq = lim u$(u) = -Di($) > according Theorem ll.181 where so is the first coefficient! 10 ! in the 

_1 

series (|2.13[) . Therefore, $(u) ~ u _1 asu-> +00. This implies the following: -F^u) ~ — u 2 asm +00. 
Thus, F2(u>i + e) > for all sufficiently small e > and -F^w) < for all sufficiently large positive u. 
Consequently, F2 has an odd number of zeroes, counting multiplicities, in the interval (w;,+oo). So F2 
has at least I zeroes in the interval (wi, +00), and it has no zeroes in (0, Wi). 

Since the first equation (|4.7[) has exactly I positive roots and the second equation (|4.7|) has at least I 
positive roots, the equation (|4.6p has at least 11 real roots, all of which are the zeroes of the polynomial p. 
But degj? = 21 by assumption, therefore, the second equation (14.71) also has exactly I positive roots. 
Moreover, it has exactly one simple root in each of the intervals (wj, £,), i = 1, ...,/ — 1, and exactly one 
simple root in the interval (w;, +00). Thus, denoting the positive roots of the second equation (|4.7p by Q , 
i = 1, . . . , I (Ci < 0), we have 

Wi < C1 2 < 6 < ^2 < Cl < 6 < ■ • • < 6-1 < W < (?. (4.9) 

Now from (gSJ|-(g3J| we obtain 

< mi < -Ci < M2 < -C2 ■ ■ • < Hi < -Cl- (4.10) 

Recall that \x% > are the positive roots of the equation (J4.6I) . and Ci < are the negative roots of 
the equation (J4.6I) . Thus, all roots of the equation (|4.6[) (and therefore, all zeroes of the polynomial p) are 
real and simple. Denote them by Xi and enumerate such that 

0<|A 1 |<|A 2 |...<|A„|. 

Then we have A2i-i = [i.% > and A 2 i = C« < 0, f = 1, . . . , I, and the inequalities f|4.10[) imply (|4. 1|) . Thus, 
p is a self-interlacing polynomial of type I. 

In the same way, one can show that if n = 11 + 1, the function $ is an i?-function with positive poles, 
and lim $(u) < 0, then the polynomial p is self-interlacing of type I. □ 

u— >±oo 

Remark 4.4. Let us note that if the even and odd parts, po and pi, of a given polynomial p have positive 
interlacing zeroes, then the polynomials po(z 2 ) and zp\(z 2 ) have real interlacing zeroes, and therefore 
the polynomial p(z) = po(z 2 ) + zpi(z 2 ) has real zeroes by Theorem 11.181 and Corollary 11.191 as a linear 
combination of polynomials with real interlacing zeroes. However, this notice does not help to investigate 
the self-interlacing property of polynomials. 

Remark 4.5. In the proof of Theorem l4.3[ we also established that the squares of both positive and negative 
zeroes of the self- interlacing polynomial p interlace both zeroes of po and zeroes of p\ . 

Theorem 14.31 allow us to use properties of i?-functions in order to obtain additional criteria of self- 
interlacing. At first, consider the Hankel matrix S($) = ||si+j||^° 7= o constructed with the coefficients of 
the series (|2.13[) . From Theorems II . 181 14.31 and 11.111 and from Corollary 11.251 we obtain the following 
self-interlacing criteria: 



10 Recall that s_i = whenever n = 11. 
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Theorem 4.6. Let a real polynomial p be defined by (|2.1[) . The following conditions are equivalent: 

1) the polynomial p is self-interlacing of type I; 

2) the following hold 

s-i < for n = 2l + l, 

r>i(*)>0, j = l,..., I, (4.11) 

Bi(*)>o, i = i,...,/, 



where I = 

3) the matrix S(^>) is strictly totally positive of rank I, and s_i < whenever n = 21 + 1. 
T7ie determinants Dj($) and Dj(Q) are defined in (|1.5[) and (11.81) . respectively. 
Proof. According to Theorem l4.3l the polynomial p is self-interlacing of type I if and only if its associated 



function $ is an i?-function with exactly I 



poles, all of which are positive, and lim $(u) = s_i sj 0. 

u— ^±00 



rm 

.2. 

By Theorem 11.181 and Corollary 11.251 this is equivalent to the inequalities (I4.11[) . But these inequalities 

are equivalent to the strictly total positivity of the matrix 5($) by Theorem II. Ill □ 



From this theorem we obtain the following criterion of self-interlacing, which is an analogue of the Hur- 
witz stability criterion. 

Theorem 4.7. Given a real polynomial p of degree n as in (j2.ll) . the following conditions are equivalent: 

1) p is self-interlacing of type I; 

2) the Hurwitz determinants Aj(p) satisfy the inequalities: 

A n _ 1 {p)>0, A n _ 3 (p) > 0, . . . , (4.12) 

rn+l] r ra+1 ] 

(-l)L 2 JA n (p)>0,(-l)L 2 i 1 A„_ 2 (p)>0,... (4.13) 

3) the determinants rjj{p) up to order n + 1 satisfy the inequalities: 

T) n (p) > 0, rin-2(p) > 0, . . . , r n (p) = a >0, (4.14) 

(-1)1 2 \ Vn+1 (p)>0,(-l)i 2 J 1 r)n _ 1 ( p ) > 0,... (4.15) 

Proof. We establish the equivalence 1) -^=^> 2) while the equivalence 2) <=>■ 3) follows from (|2.12p . 

By Theorem l4.6l p is self-interlacing of type I if and only if the inequalities (J4.11I) hold. Now from (|4. 1 1[) 
and (j2.14[) — (|2.15[) we obtain that p is self-interlacing of type I if and only if the inequalities (|4.12[) hold and 

for n = 2l (-l>»'Ajy(p) > 0, j = l,...,l; 

for n = 2l + l (-iy+ l A 2j+1 (p)>0, j = 0,l,...,l, 

that is equivalent to (|4.13p . D 

Note that the inequalities (|4.13[) and (14.15)) are equivalent to the following ones: 

n — 1 
A„_2i(j?)A„_2i_ 2 (p) < 0, and r) n -2i+i(p)r)n-2i-i(p) < 0, i = 0, 1, . . . , — — 

where A (p) = Vo{p) = a A_i(p) ee 1. 

Analogues of Theorems I3.7H3.9I will be established in Section 14.31 using a connection between Hurwitz 
stable polynomials and self-interlacing polynomials of type I. 



28 

4.2 Interrelation between Hurwitz stable and self-interlacing polynomials 

Comparing Theorems 13.31 and 14.31 one obtains the following fact, which has deep consequences and 
allows us to describe a lot of properties of self-interlacing polynomials. 

Theorem 4.8. A polynomial p{z) = po(z 2 ) + zp\(z 2 ) is self-interlacing of type I if and only if the poly- 
nomial q(z) — p(~z 2 ) — zp\(z 2 ) is Hurwitz stable. 

Proof. By Theorem l4.3[ p is self-interlacing if and only if the function $(u) = — — — - is an i?-function with 

Po(u) 

only positive poles, and lim $(u) ^ that is equivalent to the fact that the function \fr(u) = -. r is 

u->-±oo po[ — u) 

an i?-function with only negative poles, and lim ty(u) ^ 0. According to Theorem l3.3[ this is equivalent 

u— >±oo 

to the polynomial q being Hurwitz stable, as required. □ 

Remark 4.9. Thus, we have that there is one-to-one correspondence between self-interlacing and Hurwitz 
stable polynomials. Given a real Hurwitz stable polynomials, we should change sings of some its coefficients 
to obtain a self-interlacing polynomial, and vise versa. 

As an immediate consequence of Theorem 14.81 we obtain the following analogue of Stodola's theorem, 
Theorem [ 



Theorem 4.10. If the polynomial 

p(z) = a z n + a\z n ~ l + ■ ■ ■ + a n , oi, . . . , a n G R, a > 0, (4-16) 

is self-interlacing of type I, then 
for n — 11, 

(-1) 2 Oj ->0, j = 0,l,...,n; 
for n = 2^ + 1, 

(-1) 2 aj >0, j = 0,l,...,n. 
Proof. Since p is self-interlacing of type I, by Theorem 14.81 the polynomial 11 ! 

n(w+l) 

q(z) = (-1) - 5— [p (~z 2 ) - z Pl (-z 2 )} = b a z n + b^- 1 + b 2 z n - 2 + ■■■ + b n , b Q = a a > 0, (4.17) 
is Hurwitz stable. It is easy to see that 

6j = (-l) 2 oj, j = 0,l,...,n for n = 2l; (4.18) 

and 

3(3 + 1) 

bj = (-l) 2 a ,, j = 0,l,...,n for n = 2l + l. (4.19) 

By Stodola's theorem, Theorem 13. 2[ all bj are positive. D 

Thus, a necessary form of a self-interlacing polynomial with positive leading coefficient is as follows: 
for n = 21 

P (z) = b z n + b lZ n - 1 - b 2 z n - 2 - b 3 z n - 3 + b 4 z n ~ 4 + b 5 z n - 5 - b 6 Z n ~ 6 -..., b , . . . , b n > 0, 
for n = 21 + 1 

p(z) = b Q z n - hz™" 1 - b 2 z n - 2 + b 3 z n - 3 + biz" 1 ^ - b 5 z n - 5 - b 6 z n ~ 6 + ..., b ,...,b n >0. 

We also notice that if the polynomial p defined in (14.16)) is self-interlacing, then 

a„_ ia „ < 0, (4.20) 

since *(0) = -$(0) = -pi(0)/po{0) = -a„_i/a„ > (see the proof of Theorem |4~TU|) . 



n(n+l) 

11 We put the factor (—1) 2 in order to make the leading coefficient of the polynomial q equal to aa > 0. 
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Remark 4.11. In fact, Theorem 14.101 was proved in [5] by another methods. 

Let us point out at one more interesting connection between Hurwitz stable and self-interlacing polyno- 
mials. Let the polynomial p be self-interlacing. Iip(z) = po(z 2 )+zpi(z 2 ), then p(iz) — po(—z 2 )+izpi(— z 2 ), 
where i = \J — 1. Consequently, the Hurwitz stable polynomials q(z) = pq(—z 2 ) — zpi(-z) can be repre- 
sented as follows 

p(iz)+p(-iz) p(iz)-p(-iz) 1 + i . 1 - i 

q( z ) = o o- = P»«hr+K-«h-' 

2 2i 2 2 

So one can establish the following theorem. 

Theorem 4.12. Let the polynomial p{z) = po(z 2 ) + zp\(z 2 ) be self-interlacing of type I and let its dual 
polynomial q{z) = po(—z 2 ) — zpi(—z 2 ) be the Hurwitz stable polynomial associated with p. Then 

p{\) = <=*> argg(iA) = - or — ; 

and, respectively, 

q(fi) = <^=» argp(i/x) = - or — . 

Ad (A 2 ) 
Proof. In fact, p(A) = if and only if = 1. At the same time, g(zA) = po(A 2 ) — iApi(A 2 ). 

Api(A2)\ ^5^ 

— y-TT— = arctan 1 = — or — - 
Po(A 2 ) J 4 4 

The second assertion of the theorem can be proved analogously. □ 



Consequently, argg(iA) = arctan ( ) = arctan 1 = — or — 



4.3 Lienard and Chipart criterion, Stieltjes continued fractions and the signs 
of Hurwitz minors 

Now we consider the polynomials p(z) — po(z 2 ) + zpi(z ) and q(z) = po(—z 2 ) — zpi(-z) and their 

associated functions $(u) = — — — and $(u) = ; r. Let 

Po(u) Po{-u) 

So Si S-2 S3 

* «) = s -i + - + 4 + rl + ~, 



u u 2 u 3 U A 



Then we have 



Pl(—U) . . *0 . *1 , *2 , *3 *3 , S Si S 2 S 3 S3 

7 T =t-l-\ 1 2 H 3 + _ t H ^ H = _S -1 H 2 H 3 4 H 5 

po(~w) w w w° w w° w w w° w w° 



Thus, the connection between Sj and tj is as follows 

^ = (-1)^-, .7 = -1,0,1,2,... 

Consider the two infinite Hankel matrices S — ||sj+j||q° and T — ||ti+j||o°' 
In [16] there was proved the following formula. 

Theorem 4.13. Minors of the matrices S and T are connected as follows 

t( 11 h ■■ im )=(-l)£ ik+ £i jk s( i } l2 ••■ im ), m = l,2,... 

\3l ]2 ■■■ 3m J \3l J2 ■■■ 3m) 

From this theorem one can easily obtain the following consequence. 

Corollary 4.14. Let rational functions $(u) and \&(u) be such that ty(u) — — <fr(— u). Then the following 
equalities hold 

D 3 -(*) = D i (*), j = 1,2,... 

D,($) = (-lySjW, 3 = 1,2,... 
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From this corollary and from (I2.14j) - (|2.15[) one obtains 

Ti 

A n+ i_ 2i (g) = A n+ i_2j(p), i = !,•••, ^ , (4.21) 



2 



n 



A„_ 2 ,(< z ) = (-l)[ n 2 1 ] J A n _ 2j -(p), .7 = 0,1, ... , [^J , (4.22) 

where the polynomials p and g are defined in (|4.16p and (|4.17|) . respectively. 

Now we are in position to prove analogues of Theorems 13.71 1331 for self-interlacing polynomials. 

Theorem 4.15. The polynomial p defined in (|4.16p is self-interlacing of type I if and only if 

rn+l] [ n+1 ] fn±11 , 

(-l)L 2 \a n >0, (-l)L 2 J 1 o„-2>0, (-l)L 2 J 2 a „-4>0, ... (4.23) 



rn+i] r 7i+ 1 ] r n+i i 9 

(-l)L 2 J 0n >o, (-l)L 2 J 1 On _ 1 >0, (-l)L 2 J 2 a „_ 3 >0, ... (4.24) 

and one o/ i/ie following two conditions holds 

A„_i(p) > 0, A„_ 3 (p) > 0, A„_5(» > 0, . . . ; (4.25) 

r n+i i r 21+11 i rii+ii o 

(-l)L 2 JA„(p)>0, (-l)L 2 J 1 A„_ 2 (p)>0, (-l)L 2 J 2 A„_ 4 (p) > 0, . . . (4.26) 

Proof. If the polynomial p is self-interlacing of type I, then by Theorems 14.71 and 14.101 all the condi- 
tions ([4T23]) -(|4T26 )l hold. 

Let the conditions (|4.23[) and (|4.25|) hold. Consider the polynomial q defined in (|4.17|) . From (|4.23|) 
and (j4TT8l) - (J4TT9l) it follows that b n > 0, b n -2 > 0, 6„_ 4 > 0, . . ., and from (J4T251) . (f2T4|) - (f2~T5l) and (J4T2T1) we 
obtain that A„_i(q) > 0, A n _3(g) > 0, . . . Thus, the polynomial q satisfies the condition 1) of Theorem 13. 71 
Therefore, g is Hurwitz stable, so p is self- interlacing of type I according to Theorem 14.81 

Analogously, using (j4~T8| - (l4T9l) . (l2~T4l) - ([2~T5| . (|4T2T]) - ((4T22|) and Theorem [37] one can show that the 
conditions (|4~J4l) and ((4725]) . or (|4~23| and (|4~26|) . or (|4~24l) and (|4~26|) imply Hurwitz stability of the 
polynomial g, so by Theorem 14. 8[ p is self- interlacing of type I. □ 

The following theorem presents a relation between self-interlacing polynomials and continued fractions 
of Stieltjes type. 

Theorem 4.16. The polynomial p of degree n defined in (|4.16[) is self-interlacing of type I if and only if 
its associated function $ has the following Stieltjes continued fraction expansion: 

$(u) = c + , with (-l) 4 " 1 c l >0, i = l,...,2l, (4.27) 

c\u H 



CI 



1 



C 3 U 



1 



1 

c%i-\U-\ 

c 2 ; 



where cq = if n is even, and cq > if n is odd, and I as in (|2.2|) . 

Proof. In fact, by Theorem 14.31 the polynomial p is self-interlacing of type I if and only if the function $ 

can be represented in the form (14.51) . where —a = cq = if n = 21, and —a = cq = — < if n = 21 + 1. 

ai 
Now the assertion of the theorem follows from Theorem 1 1.1 81 and Corollary 1 1.391 □ 

Theorems 13.91 and 14.81 imply the following theorem. 
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Theorem 4.17. Let p be a self-interlacing polynomial of type I of degree n ^ 2 as in (|4.16p . Then all the 
polynomials 



n-2j 



n — i 


( 


n — i 


2 


2 



1 



+ j — 1 OiZ 



n—2j—i 



J = h 



1. 



aiso are self- interlacing of type I {if j is even) or of type II (if j is odd). 



Let again p be a self-interlacing polynomial of type Iof degree n and let q be its associated Hurwitz 

n(n+l) 

stable polynomial, that is, p(z) = po(z 2 ) + zp\(z 2 ) and q(z) = (—1) 2 [po(—z 2 ) — zpi(— z 2 )]. We are in 
a position to find an interrelation between Hurwitz minors of these polynomials. 

Let n = 21. By Theorem 14. 8[ the Hurwitz matrix of the polynomial q has the form 



n n ( q ) 



( a\ -a 3 a 5 -a 7 

a — a2 a4 — a 6 

a\ —as a$ 

ao — 02 04 



\o 















\ 



(-IpiaJ 



It is easy to see that the matrix 'H n (q) can be factorized as follows 

n n ( q ) = c n n n {p)E n , 

where the n x n matrices E n and C n have the forms 



E n 



(\ 
0-10 
10 
0-10 
1 



V 



/l 



L j) 






1 








.. 








-1 





.. 











-1 


.. 














1 .. 



(4.28) 



All non-principal minors of these matrices are equal to zero. So since ejj = (— l)- 7 1 and 

(J-l)0-2) 

(—1) 2 ; we have 



u 33 



Cy, 



n «2 



™ (i fc -l)(i fc -2) 



-En 



«1 «2 



(-1)* 



22 jfc-m 



(4.29) 



where 1 ^ ii < «2 < • • • < *m $J »l. Thus, the Cauchy-Binet formula together with (|4.29|) and (|4.28|) 
implies, for n = 21, 



*-<«>G I 



where 1 ^ 



™ fa-l)fa-2) m . 

\J\ J2 ■ • • Jr, 



«1 < «2 < ■ • • < V. 



< n, 



(4.30) 



J\ < ]2 < ■ ■ ■ < Jr. 
Let now n = 21 + 1. Then by Theorem 14.81 the Hurwitz matrix of the polynomial q has the form 

\ 

a ~ a 2 a 4 ~ a 6 ■ ■ ■ 

n _«. r,„ —n. n 



f-a\ 


«3 


—0,5 


a 7 





a 


-02 


f'4 


-a 6 . 








-Ol 


«3 


—as ■ 








Of) 


-02 


0,4 





I 











. (-!)[* 
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In this case, H n (q) also can be factorized: 

n n (q) = c n n n { P )(-E n ), 

where the n x n matrix C„ is as follows 

(1 
0-10000 
0-1000 
10 
10 
0-1 



L'n — 



V 



(4.31) 



All non-principal minors of the matrices C n and — E n equal zero. The principal minors of these matrices 
can be easily calculated: 



L n 



n «2 
k k 



S ifc(ifc-l) 

= (-!)£ 2 . 



E,, 



ll «2 
k k 



= (-I)*- 



(4.32) 



where 1 ^ k < k < • • • < *m ^ "• Thus, the Cauchy-Binet formula together with (|4.32[) and (|4.31|) 
implies, for n = 21 + 1, 



«n(ff) 



*i *a ... im \ =( _ 1) |/ fc( r 1) +£ i ^ ( p) fk k ... V, 

Jl J2 ■•■ Jm / \3l ]2 ■ ■ ■ Jr. 



(4.33) 



, -, / *1 < k < ■ ■ • < im ^ 

where 1 C . . ^ n. 

' Jl < J2 < • ■ • < 3m 



Since q is Hurwitz stable, from Theorem 13.51 and from the formulae (|4.30j) and (|4.33j) we obtain the 
following theorem. 

Theorem 4.18. Let p be a self-interlacing polynomial of type I of degree n and let H n (p) be its Hurwitz 
matrix defined in (|2. 10[) . Then 



for n = 21, 



forn = 21 + 1, 



(-1)*- 1 *-* «„(p) ) > 0; 

yi J2 ■•■ j^ 



™ ifc(ifc-l) ™ . /■ ■ 

(-1)*- 1 k=1 H n (p) ) > 0, 

\J1 J2 •■• J^ 



, , ., h < i 2 < ■ ■ ■ < i m s 
where 1 < . . ^ n. 

Jl <32 < ■■■ < Jm 



4.4 The second proof of the Hurwitz self-interlacing criterion 

Given a polynomial 

p(z) = a n z n + aiz n ~ 1 + ■ ■ ■ + a n , Oi,...,o n Si, a > 0, 
we consider the following rational functions 



R{z) = (-1)>(-^) and F(z)= [ 



p(z) 



p(z) 



i2(*) (-l)»p(-*) 



(4.34) 



(4.35) 
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and expand them into their Laurent series at oo: 
R(z) = 1 



£0 

z 



■2 



£2 
-3 



f(z) = i + ii + 4 + 



«a 



As we mentioned in Section ITTTl ranks of the matrices S = ||Sj+fc||o° an d T = \\tj+k\\c? are equal to the 
number of poles of the functions R and F, respectively. It is easy to see that rank of each matrix equals n 
if the polynomials p(z) and p(~z) have no common zeroes. In the rest of this section we consider only 
such polynomials, so in this section, ranks of the matrices S and T always equal n. Denoting 



q(z) = (-l)X-z), 



we obtain R 



q 



and F 



P 
P q' 

Lemma 4.19. For the functions R and F defined in (|4.35l) . the following formultE are valid: 



V 2j (p,q) = a%D 3 {R) = (-l)^Va A^ 1 (p)A ] (p), j = 1,2,. 



and 



V 2j (q,p) = a^D 3 {F) = (-l)^Va Aj-i(p)Aj(p), j = 1,2, 



(4.36) 



(4.37) 



where the determinants ^2j(p,q) are the leading principal minors of the matrix %2n+i(p, q) of order If, 
and Aj (p) are defined in (|2.1ip , Aq (p) = 1 . 



Proof. We prove the formulae (|4.36|) . The formulae (I4.37[) follow from (|4.36p since V2j(q,p) 
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Using this lemma it is easy to prove the equivalence of the conditions 1) and 2) of Theorem 14. 71 
Theorem 4.20. The polynomial p defined in (14.341) is self-interlacing of type I if and only if 

A„_i(p) > 0, A„_ 3 (p) > 0, . . . , (4.38) 



(-l)L 2 JAn(p)>0,(-l)L 2 . 

where the Hurwitz minors A^p) are defined in (|2.11[) . 



i»i-2 



(p) > o, 



(4.39) 



Proof. In the proof of Theorem 14. 3[ it was established that if p is self-interlacing of type I, then the 

7)( — Z i 

function — — is an i?-function with exactly n poles. 

p(z) 



Let n = 21 + 1. Then the function R defined in (|4.35[) is an _R-function. By Theorem 11.181 we have 
Dj(R) > for j = 1, . . . ,n. These inequalities together with the formulae (|4.36[) imply 



(-l) m i 11 A^ 1 (p)A ] (p)>0, 



J = 1, • 



.,71. 



(4.40) 



Multiplying the inequalities (|4.40[) for j — 2m and j — 2m — 1, we obtain 



A| m _ x A 2m A 



2m-2 



>0. 



Consequently, the minors A 2 j(p), i = 1, . . . , I, are positive, so the inequalities (|4.38|) are proved for odd n. 
If we multiply the inequalities (J4.40I) for j = 2m and j = 2m + 1, we get 

-AL(p)A 2m _ 1 (p)A 2m+1 ( P ) > 0. 



These inequalities imply (|4.39|) for odd 
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The converse assertion can be proved in the same way. That is, the inequalities (|4.38l) - (|4.39[) imply the 
inequalities (|4.40|) which in turn imply the inequalities Dj(R) > 0, j = 1, . . . ,n, according to (|4.36[) . By 

Theorem 11.181 the function ^— is an i?-function, so p is a self-interlacing polynomial of type I (see 



p(z) 



the proof of Theorem [4 



If n is even, the inequalities (|4.38[) - (|4.39[) can be proved analogously using the function F defined 
in (J4.35P instead of the function R. □ 



From Lemma l4.19l and from Hurwitz's stability criterion (see Theorem l3.5[) it follows one more stability 
criterion. 



Theorem 4.21. 



1) The polynomial p is Hurwitz stable if and only if for the function R defined in ((4.351) , the following 
inequalities hold 

(-l) lM ^ 1 D ] (R)>0, j = l,2,...,n. 



2) The polynomial p is Hurwitz stable if and only if for the function F defined in (|4.35[) , the following 
inequalities hold 

(-l) m ^ 1 D 3 (F)>0, j = l,2,...,n. 

Finally, we establish a relationships between the Hurwitz minors Aj(p) of a given polynomial p and 
the Hankel minors Dj(R) and Dj{F) (for the definition of the minors Dj see 1(1.80 ). To do this, consider 
the functions 



hi(z) -2a x z n - 2a. 3 z"- 2 - 2a 5 z"" 4 

zR{z) = z + -tV- = z + 



Si , S 2 



S3 



zF{z) 



p(z) 
h 2 (z) 



ciQZ n + a\z n + ■ • • + a ri 



(-l)"p(-z) 



= z + 



2a lZ n + 2a 3 z n - 2 + 2a 5 z™- 4 + . 
a z n - a x z n - x -\ h {-l) n a n 



z + s H 1 — tH — ? 

z z z* 



= z + t + — + ^ + -i 



(4.41) 



(4.42) 



where So = to = 1, and note that Dj{R) — Dj(zR) and Dj(F) = Dj(zF). This allows us to establish 
the following lemma. 



Lemma 4.22. Given the polynomial p, for the functions R and F defined in (|4.35|) . the following rela- 
tionships hold 

V 23 (p,h 1 )^a^D ] (R) = (-l) m ^ 1 VA 2 Jp), j -1,2,..., (4.43) 



Vjy(p,ft 2 ) = a^D 3 {F) = (-l)-r-2*Aj(p), j = 1,2, 



(4.44) 



where the determinants 1 ^ W 2 j (p, hi) and V 2 j (p, h 2 ) are the leading principal minors of order 2j of the ma- 
trices H 2n +i{p, hi) and H 2n +i{p,h 2 ) respectively, and Aj(p) are the Hurwitz minors of the polynomial p, 
A (p) = l. 

Proof. Assuming Qj = for j > n, from (|4.41[) we obtain 



V 2j (p,h 1 ) = {-2Y 
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2 The polynomials hi and hi are denned in I I4.41H and 114.421 1. respectively. 
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= (-l) li ^2^A J 2 ( P ). 



The relationships (I4.44[) can be proved analogously. 
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Thus, we obtain that for any real polynomial p of degree n with nonzero Hurwitz minors and, in 
particular, for Hurwitz stable and for self-interlacing polynomials, the following inequalities hold 



{-l^D^R) = (-l^D^F) > 0, j = 1,2, 



..,n. 



Last, let us note that the formulae (|4.36l) - (|4.37|) and (I4.43[) - (|4.44[) can be obtained (overcoming certain 
difficulties) from some theorems of the book [41] . But it was more simple to deduce them directly as we 
did in Lemmata 14.191 and 14.221 



4.5 Almost and quasi- self-interlacing polynomials 

In this section, we describe polynomials, which are dual (in the sense of Theorem l4.8p to the quasi-stable 
polynomials. Because of the mentioned duality we just give the definition of these polynomials. 

Definition 4.23. A polynomial p(z) — po(z 2 ) + zpi(z 2 ) of degree n is called quasi-self-interlacing of 
type I (of type II) with degeneracy index m, ^ m ^ n, if the polynomial po(— z 2 ) — zpi{—z 2 ) (rcsp. the 
polynomial Pq(— z 2 ) + zpi(—z 2 )) is quasi-stable with degeneracy index m. 

In other words, the polynomial p is quasi-self-interlacing if non-common zeroes f the polynomials p(z) 
and p{— z) are real, simple and interlacing. 

Obviously, all results regarding quasi-stable polynomials can be easily reformulated for the quasi-stable 
polynomials. So we leave such reformulation to the reader. 

In Section [S] we use the quasi-stable polynomials with degeneracy index 1, that is, the polynomials of 
the form p(z) = zq(z), where q(z) is a self-interlacing polynomial. Such polynomials we call almost self- 
interlacing polynomials. Note that if p(z) = po(z 2 ) + zpi(z 2 ) is almost self-interlacing, then the function 

$ = — is an i?-function with all poles positive except one, which is zero. 
Po 
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4.6 "Strange" polynomials 

This section is devoted to short description of some numeric results on the distribution of zeroes of 
some polynomials closely connected to the Hurwitz stable and self-interlacing polynomials. We call these 
polynomials " strange" because of their curious zero location. 

Let p(z) — Po(z 2 ) + zpi{z 2 ) be a Hurwitz stable polynomial of degree n. Consider the polynomial 
q(z) = po(—z 2 ) + zpi(z 2 ). It is evident that q(0) = p(0) ^ 0. Calculations with Maple Software showed 



that the polynomial q has exactly 



n + 1 



simple zeroes in the open right half-plane and exactly 

simple zeroes in the open left half-plane, but it has no zeroes on the imaginary axis. At least one zero of 
q is nonreal. 

Let us denote by /_ti , /12 , ■ • ■ the distinct absolute values of the zeroes of the polynomial q lying in 
the open left half-plane such that 

< Mi < M2 < Ms < • • • 

And denote by Ai, A2, • . • the distinct absolute values of the zeroes of the polynomial q lying in the open 
right half-plane such that 

< Ai < A 2 < A 3 < . . . 

So numerical experiments showed that 

< Ai < m < A 2 < p 2 < A3 < • ■ ■ 

The polynomial po(z 2 ) + zpi(—z 2 ) possesses similar properties according to calculations. 

However, it is not known for sure if all the polynomials of the type po(— z 2 ) + zpi{z 2 ), where the 
polynomial pq(z 2 ) + zpi(z 2 ) is Hurwitz stable, have the zero location described above. 

4.7 Matrices with self-interlacing spectrum 

In this section, we consider some classes of real matrices with self-interlacing spectrum and develop a 
method of constructing such kind of matrices from a given totally positive matrix. But at first, we recall 
some definitions and statements from the book [ID] . 

Definition 4.24 (p~PJ). A square matrix A = ||a,j||" is called sign definite of class n if for any k ^ n, all 
the non-zero minors of order k have the same sign Sk- The sequence {ei, £2, • ■ • , £ n } is called the signature 
sequence of the matrix A. 

A sign definite matrix of class n is called strictly sign definite of class n if all its minors are different 
from zero. 

Definition 4.25 ([ID]). A square sign definite matrix A = || Oy ||" of class n is called the matrix of class n + 
if some its power is a strictly sign definite matrix of class n. 

Note that a sign definite (strictly sign definite) matrix of class n with the signature sequence e± = 
e-2 = ■ ■ ■ = £ n = 1 is totally nonnegative (strictly totally positive). Also a sign definite matrix of class n + 
with the signature sequence £\ = £2 = ■ ■ ■ = s n — 1 is an oscillating matrix (see [10]). that is, a totally 
nonnegative matrix whose certain power is strictly totally positive. It is clear from the Binet-Cauchy 
formula that the square of a sign definite matrix is totally nonnegative. 

In [TO] it was established the following lemma. 

Lemma 4.26. Let the matrix A = ||ay||™ be totally nonnegative. Then the matrices B — ||a n — i+i,j||™ and 
C = ||ai,n-j+i||i are sign definite of class n. Moreover, the signature sequence of the matrices B and C is 

as follows: 

fc(fc-i) 
£ fe = (-l) 2 , fc = l,2,...,n. (4.45) 

Note that the matrices B and C can be represented as follows 

B = J A, and C = AJ, 
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where 



J = 



/O . 
. 


. 

. 1 


1\ 




1 . 
\l . 


. 
. 






(4.46) 



It is easy to see that the matrix J is sign definite of class n (but not of class n + ) with the signature 
sequence of the form (|4.45[) . So by the Binet-Cauchy formula [TT] we obtain the following statement. 

Theorem 4.27. The matrix A = ||ay||™ is totally nonnegative if and only if the matrix J A (or the 
matrix AJ ) is sign definite of class n with the signature sequence (J4.45P . The matrix J is defined in (J4.46I) . 

Obviously, the converse statement is also true. 

Theorem 4.28. The matrix A = 



Ojjll" is a sign definite matrix of class n with the signature se- 
quence (|4.45[) if and only if the matrix J A (or the matrix AJ) is totally nonnegative. 

In the sequel, we need the following two theorems established in the book [TO] . 

Theorem 4.29. Let the matrix A = ||o,j ;||" be a sign definite of class n + with the signature sequence £&, 
k = 1, 2, . . . , n. Then all the eigenvalues A&, k = 1, 2, . . . ,n, of the matrix A are nonzero real and simple, 
and if 

|Ai|>|A 2 |>...>|A„|>0, (4.47) 



then 



sign X k 



£fc 

Eft— 1 



k = 1,2, ...,n, £ = 1. 



(4.48) 



Theorem 4.30. A totally nonnegative matrix A = ||ay||™ is oscillating if and only if A is nonsingular 
and the following inequalities hold 



l 3,j + 



i>0, 



and 



Oj+ij > 0, j = 1,2,. 



Let us introduce the following definition. 

Definition 4.31. A matrix A is said to have the self-interlacing spectrum if its eigenvalues are real and 
simple and satisfy the following inequalities 



or 



Ai > -A 2 > A 3 > . . . > (-l) n_1 A„ > 0, 

- Ai > A 2 > -A 3 > . . . > (-1)"A„ > 0. 
Now we are in a position to complement Lemma 14.261 



(4.49) 



(4.50) 



Theorem 4.32. Let all the entries of a nonsingular matrix A — Ha^H" be nonnegativq^j and for each i, 
i = 1,2, . . . , n — 1, there exist number n and r%, 1 ^ r\, r-i ^ n, such that 



&n—i.r\ ' On+1 — r*i,i ^ U, Gn+1 — i,r2 ' ^n+ 



1— r 2 ,i+ 



i >o 



[Or ^i,n+l — r± ' &ri,n—i -^ ") ^i+l,n+l— r% ' ^T2, 



n+1— i 



>o). 



(4.51) 
(4.52) 



The matrix A is totally nonnegative if and only if the matrix B = J A = ||&;j||i (or, respectively, the matrix 
C = AJ) is sign definite of class n + with the signature sequence defined in (|4.45j) . Moreover, the matrix B 
(or, respectively, the matrix C) possesses a self-interlacing spectrum of the form (J4.49I) . 



'The matrices with nonnegative entries are usually called nonnegative matrices. 
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Proof. We prove the theorem in the case when the condition (14.51)) holds. The case of the condition (|4.52[) 
can be established analogously. 

Let A be a nonsingular totally nonnegative matrix and the condition (I4.51|) holds. From Theorem l4.27l 
it follows that the matrix B = J A is sign definite of class n with the signature sequence (|4.45p . In order to 
the matrix be be sign definite of class n + it is necessary and sufficient that a certain power of this matrix 
B be strictly sign definite of class n. Since the entries of the matrix J have the form 

, n _) 1, i = n+l-j; 

{ 0, i^n + l-j; 

the entries of the matrix B can be represented as follows: 

n 

Oij = 2_^(J)ik a kj = O-n+l-iJ- 
fe=l 

Consider the totally nonnegative matrix B 2 . Its entries have the form 

n n 

(B )ij = 2_^bikbkj = > Qn+l-i,fcQra+l-fc,j- 

fc=l fc=l 

From these formulas and from (|4.51[) it follows that all the entries of the matrix B 2 above and under the 
main diagonal are positive, that is, (B 2 )i,i + i > (B 2 )i + i : i > 0, i — 1,2, ...,n— 1. By Theorem 14.301 
B 2 is an oscillating matrix. According to the definition of oscillating matrices, a certain power of B 2 is 
strictly totally positive. Thus, we proved that a certain power of the matrix B is strictly sign definite, so 
B is a sign definite matrix of class n + with the signature sequence (J4.45I) according to Lemma 14.261 By 
Theorem 14.291 all eigenvalues of the matrix B are nonzero real and simple. Moreover, if we enumerate the 
eigenvalues in order of decreasing absolute values as in (|4.47p , then from (|4.48|) and (I4.45[) we obtain that 
the spectrum of B is of the form (|4.49|) . 

The converse assertion of the theorem follows from Theorem 14.271 □ 

Remark 4.33. For n = 11 + 1 (for n = 21), the characteristic polynomial of the matrix B described in 
Theorem 14.321 is self-interlacing of type I (of type II). 

Remark 4.34. If the matrix —A is totally nonnegative and the conditions (|4.51[) (or the conditions fl4.52[) ^ 
hold, then the matrix B = J A (or, respectively, C = AJ) has a spectrum of the form (|4.50|) . 

Remark 4.35. One can obtain another types of totally nonnegative matrices which result matrices with self- 
interlacing spectra after multiplication by the matrix J. To do this we need to change the conditions (|4.51|) - 
(|4.52[) by another ones such that, for instance, the matrix B A , or B 6 , (or B 8 etc.) becomes oscillating. 

Theorem l4. 321 implies the following corollary. 

Corollary 4.36. A nonsingular matrix A with positive entries is totally nonnegative if and only if the 
matrix B = J A (or the matrix C = AJ) is sign definite of class n + with the signature sequence (|4.45[) . 
Moreover, the matrix B has a self-interlacing spectrum of the form (|4.49|) . 

Consider a particular case of conditions (14.5ip . Suppose that all diagonal entries if the matrix A be 
positive ctjj > 0, j = 1, 2, . . . , n. It is easy to see that in this case the conditions (|4.51|) hold if <Xj,j+i > 0, 
j = 1, 2, . . . , n — 1. If all remaining entries of the matrix A are nonnegative, then by Theorem 14.321 the 
matrices B = J A and C = AJ have self-interlacing spectra if A is totally nonnegative. 

If all other entries of the matrix A (that is, all entries except ajj and a JJ+ i, which are positive) equal 
zero, then A is a bidiagonal matrix with positive entries on and under the main diagonal. Clearly, A is 
totally nonnegative. Then by Theorem 14.321 the matrix B = J A has a self-interlacing spectrum. Thus we 
obtain the following statement established in [15) : 
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Theorem 4.37. Any anti-bidiagonal n x n matrix with positive entries 



B = 



f° 
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b n -l 
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o / 



(4.53) 



has a self-interlacing spectrum of the form (I4.49|) . Here the entries bj > 0, j = 2, 3, . . . , n, lie above the 
main diagonal, the entries Cj > 0, j = 2, 3, . . . , n, lie under the main diagonal and the only entry on the 
main diagonal is a\ > 0. 



K = 



In [T5] it was also proved that the spectrum of the matrix (|4.53p coincides with the spectrum of the 
following tridiagonal matrix 

/an 6 2 ... 0\ 



b n 

Cn 0/ 

It is well-known 10 that the spectrum of this matrix does not depend on the entries bj and Cj separately. It 
depends on products bjCj, j = 2, 3, . . . ,n. So in order to the matrices (|4.53[) and K to have self-interlacing 
spectra, it is sufficient that the inequalities a\ > and bjCj > 0, j = 2, 3, . . . ,n, hold. 
Finally, consider a tridiagonal matrix 



C2 





b 3 





C3 














V° 









/oi 


6i 


. 





\ 


Cl 


«2 


&2 • 











C2 


a 3 • 














. 


ln-1 


b n -l 


Vo 





. 


• Cn-1 


a n } 



Mj = 



where a,k, &fe, c^ S M and Ckbk ^ 0. In [TU], there was proved the following fact. 

Theorem 4.38. The matrix Mj is oscillating if and only if all the entries bk and Ck are positive and all 
the leading principal minors of Mj are also positive: 



a x > 0, 
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(4.54) 



This theorem together with Theorem 14.321 implies the following statement. 
Theorem 4.39. The anti-tridiagonal matrix 
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where a,j,bj,Cj > for j = l,2,...,n — 1, is sign definite of class n + and has a self-interlacing spec- 
trum (|4.49p i/ and only if the following inequalities hold: 

(-V^Ajf ,] , J. , ••■ fc )>0, k=l,2,...,n. (4.55) 

\n + 1 — k n + 2 — k . . . nj 

Proof. If the matrix Aj is sign definite of class n + and has a self-interlacing spectrum of the form (14.491) . 
then according to Theorem 14.291 the signs of its nonzero minors can be calculated by the formula (|4.45l) . 
This implies the inequalities (|4.55|) . 

Conversely, let the inequalities (|4.55[) hold. Then we have that the inequalities (J4.54I) hold for the matrix 
Mj = JAj. By Theorem 14.381 the matrix Mj is oscillating and, in particular, totally nonnegative. 
Now notice that (Mj)u > 0, i — l,...,n, and (Mj)k,k+\ > 0, k = 1, . . . , k — 1, so Mj satisfies the 
condition (|4.51[) of Theorem 14.321 Therefore, the matrix A j is sign definite of class n + and has a self- 
interlacing spectrum of the form (J4.49I) . □ 

5 Generalized Hurwitz polynomials 

In this section, we describe the class of real polynomials whose associated function $ is an i?-function. 
In particular, this class includes all Hurwitz stable and self-interlacing polynomials as extremal cases. 

5.1 General theory 

Let us consider a polynomial 

p(z) = a z n + axz n ~ + ■ ■ ■ + a n , Oi,...,o„el, a > 0. (5.1) 

Definition 5.1. The polynomial p is called generalized Hurwitz polynomial of type I of order k, where 

In + ll 
1 ^ k ^ — - — , if it has exactly k zeroes in the closed right half-plane, all of which are nonnegative and 

simple: 

sj m < \i 2 < ■■■ < Vk, (5.2) 

such that p(— Hi) ^ 0, i = 1, . . . , k, and p has an odd number of zeroes, counting multiplicities, on each 
interval (— /Ufc,— /^fc-i), • • • , (— M3, — 1^2)1 (~ l J >2i~IJ'i)- Moreover, the number of zeroes of p on the inter- 
val (— fii,0) (if any) is even, counting multiplicities. The other real zeroes lie on the interval (—00, — /Xfc): 
an odd number of zeroes, counting multiplicities, when n — 21, and an even number of zeroes, counting 
multiplicities, when n = 21 + 1. All nonreal zeroes of p (if any) are located in the open left half-plane of 
the complex plane. 

Definition 5.2. If p(z) is a generalized Hurwitz polynomials of type I, then the polynomial p(—z) is called 
generalized Hurwitz polynomial of type II. 

It clear that all results obtained for the generalized Hurwitz polynomials of type I can be easily refor- 
mulated for the generalized Hurwitz polynomials of type I. Thus, in the rest of the section we consider 
only generalized Hurwitz polynomials of type I unless explicitly stipulated otherwise. 

Since the generalized Hurwitz polynomials of order k have exactly k zeroes in the closed right half- 
plane, the generalized Hurwitz polynomials of order have no zeroes in the closed right half-plane, so they 

are Hurwitz stable. 

I" n _|_ \ ' 

Analogously, if p is a generalized Hurwitz polynomial of order k — without a root at zero, then 

p is self-interlacing of type I. In fact, let < [i\ < ^2 < ■ ■ ■ < Hk be its positive zeroes. Then by definition, 
p has an odd number of zeroes, counting multiplicities, (at least one) on each interval (— (ij+i, — fij), 
j = 1, . . . , k — 1. Moreover, if n — 21, then p has at least one zero on the interval (-co, —fik)- It is easy 
to see that p can not have nonreal zeroes and has exactly one simple zero on each interval (— /i J+1 , —/J,j), 
j = 1, . . . , k — 1, and on the interval (—00, —(J-k) for n — 21. Besides, the zeroes of p are distributed as 
in ([4~Tj) . 
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Thus, Hurwitz stable and self- interlacing polynomials are generalized Hurwitz polynomials of minimal 
and maximal orders, respectively. We note that generalized Hurwitz polynomials of maximal order with a 
root at zero are almost self-interlacing. 

Now we are in a position to establish the main theorem of the theory of generalized Hurwitz polynomials. 
This theorem is a direct generalization of Theorems 13.31 and 14.31 



(for 



Theorem 5.3. Let p be a given real polynomial of degree n Js 1 as in (|5.ip . The polynomial p is generalized 
Hurwitz if and only if its associated function <& defined in (|2.7[) is an R-function with exactly I = — 
ai ^ 0) or I — 1 (for ai = 0) po/ei 14 !. Moreover, 

• if n — 21 or if n = 21 + 1 with 

aoai > 0, 

then the number of nonnegative poles of the function $ equals order of p; 

• if n = 21 + 1 with 

a ai < 0, 

then the number of nonnegative poles of the function $ equals order of p minus one. 

~n+ 1" 



Proof. As we said above, for k = and k 



this theorem is true. 



Sufficiency. Let the function <I> be an i?-function with exactly r nonpositive poles, $J r ^ 

At first, assume that a\ ^ and the function $ has no pole at zero. Then $ has exactly / 
and can be represented in the following form by Theorem 1 1.1 81 



(5.3) 



(5.4) 



poles 



$(u) = a +y^ ■ 



tv, 



-+ E - Jf2 ~. (5-5) 

j — 1 J j — r+1 J 

where a £ K (a = if and only if n = 21), otj > 0,u>j > for j = 1, . . . ,1, all poles are distinct. Let the 
positive poles are enumerated as follows < Wi < L02 < ■ ■ ■ < u r . 

As in the proofs of Theorems 13.31 and 14.31 we note that the set of zeroes of the polynomial p coincides 
with the set of roots of the following equation 



$(z 2 ) = — . (5.6) 

z 

So if we suppose that p(X) — and Im A =/= (that is, A is a nonreal zero of p), then from (|5.6p we obtain 

Im$(A 2 ) 1 ImA 1 1 



Im(A 2 ) |A| 2 ImA 2 2|A| 2 ReA 



<0, 



(5.7) 



since $ is an i?-function by assumption. From (|5.7p it follows that all nonreal zeroes of p lie in the open 
left half-plane. Thus, all zeroes of p are located in the open left half-plane or on the nonnegative real 
half-axis. 

As in the proof of Theorem 14. 3[ we investigate the distribution of the real zeroes of the polynomial p. 
By assumption, |$(0)| < 00, so the polynomial p does not vanish at zero. In fact, if p(Q) — 0, then 
Po(0) = 0, but it contradicts with |$(0)| < 00. Thus, we suppose that z£l\ {0} and change the variables 
as follows z 2 = u > 0. Then for real nonzero z, the equation (15.61) is equivalent to the following system 




u >0. 



(5.8) 



We also note (as in the proof of Theorem 14. 3 1 that all positive roots (if any) of the first equation (|5.8[) are 
squares of positive roots of the equation (|5.6p , and all positive roots (if any) of the second equation 



14 Let us note that if n = 21 and <E> is an R-function, then a\ 7^ 0, so <t> has exactly I poles. 
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are squares of negative roots of the equation (|5.6p . The equation (|5.6p may have nonreal roots, but as we 
showed above all these roots are located in the open left half-plane. 

Let n = 21. Then lim Q(u) = a = 0. The set of zeroes the function F\(u) — $(u) H = coincides 

u->±oo yju 

with the set of roots of the first equation (|5.8j) . In the same way as in the proof of Theorem 14.31 one 
can show that the function F\ has exactly one simple zero, say [i\ {jx% > 0), on each interval (u>i—i,u>i), 
i = 2, . . . , r and exactly one simple zero \i\ {p,\ > 0) on the interval (0, Wi). Moreover, denoting by £i a 
unique zero 15 ! of the function $ on the interval (wj, Wj+i), i = 1, . . . , r — 1, we have 

Mi < Wi < & < /J? +1 < Wi+i, i = l,...,r-l, (5.9) 

since i*i(£j) = — p= > 0, and F\ is decreasing between its poles on the positive half-axis as a sum of 

decreasing functions. 

Thus, the first equation (J5.8I) has exactly r positive roots /zf all of which are simple and satisfy the 
inequalities (|5.9p . 

Now consider the function i^fu) = <&(-u) = whose set of zeroes coincides with the set of roots of the 

second equation (|5.8|) . In the same way as in the proof of Theorem l4.3[ one can show that the function F2 
has an odd number of zeroes, counting multiplicities, on each interval (wj, Mi+i)> i = 1, . . . ,r — 1, and on 
the interval (w r , +00). Moreover, F2(u) < on the intervals [m?+i> Wj+i), i = 1, . . . , r — 1. Let us now turn 

our attention to the interval (0,Wi). Since |3>(0)| < by assumption and -= —} —00 as u — >• +0, we 

V M 
2 ,2 

have i*2(u) — > —00 as u — >• +0. Besides, F2U1) = -Fi(u) — , so i^fwf) = < 0. Therefore, F2 has an 

V" Mi 

even number of zeroes, counting multiplicities, on the interval (0,^f). However, since F\(u) is monotone 
decreasing on (0,wi) as a sum of two monotone decreasing functions, we have F\{u) < on the interval 
(Mi,wi). Consequently, F^v) is also negative for u € (/xf,wi), so its zeroes on the interval (0,wi) (if any) 
lie, indeed, on the interval (0,^). 

Thus, we obtain that in the closed right half-plane, the polynomial p has only r zeroes, all of which are 
positive and simple: 

< [il < ^2 < • ■ ■ < Mr- 

Moreover, p has an even number of zeroes, counting multiplicities, on the interval (— Mi)0) and an odd 
number of zeroes, counting multiplicities, on each interval (— Mi+ij ~M«)> i = 1, ■ ■ . ,r — 1. Also p(— Mi) 7^ 
for i = 1, . ..,r. As we showed above, all nonreal zeroes of p are in the open left-half-plane, so p is a 
generalized Hurwitz polynomial of order r. 

Let n. = 21 + 1, and the condition (|5.3[) holds. The only difference between this case and the case n = 21 

is behaviour of the function <I> at infinity. Thus, as in the case n = 21, i<\ has a unique simple zero, say /if 

(Mi > 0), i = l,...,r, on each interval (0,u;i), (wi,W2), ..., (w r _i,u; r ). These zeroes are distributed as 

in (|5.9p . Furthermore, we have Fi(u) — > +00 as u \ w r and Fi(u) — > a = — > as u — > +00, so 

ai ■ 
Fi(u) > for u G (w r , +00) (recall that $ is a monotone decreasing function on the interval (w r , +00) and 

$(u) — > +00 asu\ tu r ). As in the case n = 21, the function i^Cw) has an even number of zeroes, counting 

multiplicities, on the interval (0, /if) and an odd number of zeroes, counting multiplicities, on each interval 

(u>j,m!+i) f° r i = l,...,r — 1. Moreover, F2 has no zeroes on the intervals (/if,wi) and (^ +1 ,u;i+i), 

i = 1, . . . ,r — 1. At the same time, we have F<z{u) — > +00 as u \ ui r and F2(u) — > a > as u — > +00. 

Since the function ^(u) is not monotone, it has an even number of zeroes, counting multiplicities, on 

the interval (w r , +00). Recall that these zeroes are squares of positive zeroes of the polynomial p on the 

interval (w r , +00). Since the other zeroes of the functions F\ and Fi are distributed as in the case n = 21, 

we obtain that p is a generalized Hurwitz polynomial of order r. 

Let now n = 21 + 1, and let a^ai < 0. This case also differs from the case n — 21 by behaviour of the 
function $ at infinity: lim $(w) = a = — < 0. As in the case n = 21, the function F\(v) — $(w) H -= 

has a unique simple zero, say /j,f (/Xj > 0), i = 1, . . . , r, on each interval (0, Wi), (wi, W2)) • • • , (w r -i, w r ). 



15 Recall that $ is an /J-function by assumption. Consequently, it has exactly one simple zero on each interval (u)j, u;;_|_i). 
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These zeroes are distributed as in (|5.9p . Furthermore, we have Fi(u) — >• +00 as u \ w r and Fi(u) — > a < 
as u-> +00. Since Fi(u) is decreasing on the interval (uv, +00) as a sum of two decreasing functions on 
this interval, it has a unique simple zero, say /^ +1 (/V+i > 0) on (u> r , +00) such that 

OJr < & < Mr+1; 

where £ r is a zerq 16 l of <& on the interval (w r ,+oo). As in the case n — 21, the function i*2(u) has 
an even number of zeroes, counting multiplicities, on the interval (0, /uf ) and an odd number of zeroes, 
counting multiplicities, on each interval (uii,/j,f +1 ) for i — l,...,r — 1. Moreover, F? has no zeroes on 
the intervals (/if, Wj), i — 1, . . . , r. Now consider the interval (u r , +00). Since F2(u) — > +00 as u \ uj r 

and -F^Mr+i) = < 0, we obtain that F2 has an odd number of zeroes, counting multiplicities, 

Mr+l 

2 

on the interval (w r , M?+i) • Moreover, i^w) = F\{u) = < for (/i,^ +1 ,+oo), because F\ is negative 

V u 
on that interval. Since the positive roots of the first equation (|5.8|l are squares of the positive zeroes 
of the polynomial p and the roots of the second equation (J5.8I) are squares of the negative zeroes of the 
polynomial p, we obtain that that p has only r + 1 zeroes in the closed right half-plane, all of which are 
positive and simple: 

< /ii < [12 < ■ ■ ■ < fJ-r+l- 

Moreover, p has an even number of zeroes, counting multiplicities, on the interval (— /Ui,0) and an odd 
number of zeroes, counting multiplicities, on each interval (— Hi+i, —fii), i = I, . . . ,r. All nonreal zeroes 
of p are in the open left half-plane, and p(—fn) =/= for i = 1, . . . , r + 1, so p is a generalized Hurwitz 
polynomial of order r + 1 . 

If n = 21 + 1, and ai = 0, then by Theorem 1 1.181 the function $ has exactly I — 1 poles and can be 
represented in the form 

r l-l 

$(u) = (3u + 1 + J2^^+ E -r-> ( 5 - 10 ) 

* — ' M — UJj A — ' U + USj 

j — 1 J ^— r+1 J 

where /) < 0, 7 £ 1, Oj > 0, w^ > for j — 1, . . . ,1, all poles are distinct. In the same way as in the case 
n = 21 + 1 with aQdi < 0, one can show that p is a generalized Hurwitz polynomial of order r + 1, since 
$(u) is also negative for sufficiently large u. 

Suppose now that the function $ has a pole at zero. Then po(0) = 0, so p(0) = a n — 0. At the 
same time, pi(0) = a„_i = p'(0) 7^ 0, since $ is an i?-function by assumption, so it has interlacing poles 
and zeroes (see Theorem 1 1.18(1 . Thus, p has a simple root at zero. Note that the set of zeroes of the 
polynomial p, but the root at zero, coincides with the set of roots of the equation (J5.6I) . 

If n = 21 or n = 21 + 1 with (|5.3[) , then as above, one can show that in the closed right half-plane p has 
exactly r zeroes, say /Xj, i = 1 . . . , r, all of which are nonnegative and simple: 

= Hi < ^2 < i*2 < • ■ • < Mr, 

and p{~Hi) 7^ for i = 2, . . . , r. Also p has an odd number of zeroes, counting multiplicities, on each 
interval (— [Ai, —fJ-i-i) an d an odd (even) number of zeroes on the interval (—00, —fir) ii n = 21 (n = 21 + 1). 
All nonreal zeroes of p lie in the open left half-plane. Thus, p is generalized Hurwitz of order r. 

If n = 21 + 1 and the condition (|5.4|l holds, then in the same way as above, one can show that p is 
generalized Hurwitz of order r + 1. 

Necessity. Let p be a generalized Hurwitz polynomial of order k. At first, suppose that p(0) =/= 0. 

Let n = 21. Then by definition, p{z) has only k zeroes in the closed right half-plane, all of which are 
positive and simple: < (Xi < [j,2 < • • • < Mfc> and p(—fn) 7^ 0- Moreover, p has an odd number of 
zeroes, counting multiplicities, on each interval (—00, —fik), (— A*fe, — Mfe-i)r • • > {~^2i — A*i)j an d it has an 
even number of zeroes, counting multiplicities, on the interval (— /xi,0). All nonreal zeroes of p lie in the 
open left half-plane. Thus, we can represent the polynomial pas a product p(z) = h(z)g(z), where h is a 



^Recall that this zero is unique, since <& is a monotone function on (cj r ,+oo). 
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self-interlacing polynomial 17 !, of degree 2k whose positive zeroes are /Uj, i = l,...,k, and g is a Hurwitz 
stable polynomial of degree 2(1 — k), which has an even number of zeroes, counting multiplicities, on each 
interval (-oo,-/i fc ), (-fa,~Hk-i),- • • , (-fa,- Hi), (~fa,0). Let 

h(z) = h (z 2 ) + zh 1 (z 2 ), 

g(z) = g Q (z 2 ) + zgi(z 2 ). 

Then by Theorem I4.3| the function H := - — is an i?-function with exactly k poles, all of which are positive. 

h 

Analogously, by Theorem 13.31 the function G :— — is an R function with exactly I — k poles, all of which 

90 
are negative. If p(z) — Po(z 2 ) + zpi(z 2 ), then 

Pi(u) = h (u)gi(u) + hi(u)g (u), (5.11) 

Po(w) = h (u)go(u) + uhi(u)gi(u). (5-12) 

Thus, 

pi(u) _ h (u)gi(u) + hi(u)g (u) _ H(u) + G(u) 



*(«) 



Po(u) ho(u)go(u) + uh\(u)g\(u) 1 + uH(u)G(u) 

Note that the function 

h i(u) 9i(u) = pi(u) 

h (u) g (u) ho(u)go(u) 

is an i?-function as a sum of i?-functions. Consequently, all zeroes of p\ are real and simple and interlace 
zeroes of the polynomial /io3o ; which has k positive simple zeroes and / — k negative simple zeroes. 
Furthermore, the function 

h (u) ugi(u) 1 n , , Po(u) , . 

fti(u) .go(w) -ff(w) hi(u)g (u) 

is an iE-function by Theorem ll.23[ since the functions iJ and G are i?-functions. Consequently, all zeroes 
of the polynomial po are real and simple and interlace zeroes of the polynomial higo (see Theorem II . 18[) . 

Thus, from (|5.13[) (I5.14[) it follows that all poles and zeroes of the function $ = — are real and simple. 

Po 
Therefore, $ can be represented as follows 

, > Pi (u) \-^ a, 
po(u) £ju-wj 



where Wj € R, Wj 7^ cjj for i =/= j, and 

Pi(^i) 
Q J = "77 — \ 



7^0, j = l,2,...,L 



According to Theorem I1.18[ now it is sufficient to establish positivity of the numbers ay to prove that $ 
is an i?-function. 

Let uij be a zero of the polynomial po- Since p(0) 7^ by assumption, we have po(0) 7^ 0, so ojj 7^ 0, 
j = 1,2, ... , L It is clear that h\(uij)go(ujj) 7^ 0, because the zeroes of the polynomial po interlace the 
zeroes of hig$. From (|5.12j) it follows that 



h (ujj) 



Uj h i(uj)9i(uj) 

o> — t \ • 

9o(uj) 



so 

/ \ 1 / s / s W 7^l( W 7")5l ( W 7") 

Pi.(uii) = hr^Muj) - 3 Jyiy 3 ' = 

9o(uj) 

= ¥S\ W u *> ~ ">&"$ = h ij^9(^)9(-^)- 



17 More exactly, h is self-interlacing of type I. 
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According to (|5.14[) . we have Po{u) = —hi(u)g (u)F(u), and F(uij) = 0, therefore 

Po(uj) = -hiiu^goiu^F'fa). 
So we obtain _ 



p' {ujj) g (uj) -hi{uj)go{uj)F'{u)j) 9a(wj)F'(wj) 

Since F is an i?-function, we have F'(u) < 0, so 9q(wj)F'(wj) < 0. Thus, we obtain that 

sign(ay) = sign (g(y/uj)g(-y/Wjj) . 

If uij < 0, then g(y/UJj)g(—y/uJj) = <?o( w j) — 0J j9i( UJ j) > 0; since go and g\ are real polynomials. So 
cij > in this case. 

Let now ujj > 0. The polynomial g is Hurwitz stable with positive leading coefficient by construction, 
so g(z) > for all z ^ 0. Therefore, sign(aj) = sign (<;(— y/Uj))- Since degpo = U deggo = I — k, 
deg/ii = k — 1, the function F defined in (15.14ft has exactly I zeroes and I — 1 poles. Consequently, the 
maximal positive zero of the polynomial po is greater than the maximal positive zero of the polynomial hi 
(recall that all zeroes of hi are positive). Besides, F is decreasing between its poles, so by (J5.12J) . we obtain 

F(0) = — - — — — — = — - — — > 0, since h is self-interlacing by construction, and (14.201) holds for it. This 
hi(0)go{0) hi(0) 

means that the function F (and the polynomial po) has a positive zero, which is less than the minimal 

positive zero of hi- Thus, the polynomial po has exactly k positive zeroes and, respectively, exactly I — k 

negative zeroes. 

Let us enumerate the positive zeroes of po as follows 

< Ull < U>2 < ■ . . < UJk ■ 

Then denoting the zeroes of the polynomial hi by jj (j = 1, 2, . . . , k — 1), we obtain 

< Wl < 71 < L0 2 < 72 < • ■ • < Wfe < 7fe_l < LO k . 

Let us now show that wi > /if. Indeed, since /ii is the minimal positive zero of h and h(0) ^ 0, 
the polynomial h does not change its sign on the interval [0,/ii). Moreover, both polynomials ho and hi 
have no zeroes on the interval [0,/xf] (see the proof of Theorem 14.3ft - Consequently, we get sign(/i(0)) = 
sign(/io(0)) = — sign(/ii(0)) = sign (/io(a*i)) = —sign (/ii(Mi)) that follows from the self-interlacing of h. 
Note that sign (po(0)) = sign (ho(0)go(0)) = sign (ho(0)), since g is Hurwitz stable with the positive leading 
coefficient, so <?o(0) > 0. The number fix is a zero of the polynomial h, therefore ho{^\) = — /ii/ii(/zf), so 
from (|5.12ft we have 

Po(Mi) = ho{ni)g (l4) ~ Mi^o(Mi).9i(^i) = MMiM-Mi)- 

By construction, the polynomial g has an even number of zeroes, counting multiplicities, on each interval 
(-oo,-/i fc ),(-/iA : ,-p fe _i),...,(-/i2,-Mi),(-A i i:0) I therefore sign(p(0)) = sign (g(-fi 3 )), j = 1,2, . . . , k, 
that is, g(—fj,j) > 0. This implies that sign (po (/if)) = sign (ho ((if )) = sign(/io(0)) = sign(po(0)). Conse- 
quently, the polynomial po has an even number of zeroes, counting multiplicities, on the interval (0,/if). 
But the positive zeroes of the polynomials po and hi interlace, and hi has no zeroes in the interval (0,/if), 
therefore po has also no zeroes on that interval, that is, u\ > (i\ ■ 

Note that (see the proof of Theorem [ 



< Mi < 71 < nt < 72 < • • • < Mfc-i < 7fc~i < A- 



Thus, we already established that < fif < u>i < 71 < fj%. Now we prove that the polynomial po 
has no zeroes on each interval (7^, /t| +1 ), j — 1, 2, . . . , k — 1. In fact, since sign (po(0)) — sign (h (0)) and 
both polynomials po and /io change their signs on each interval (77, 7j+i), J = 1, 2, . . . , fc — 1, we obtain 
sign (po(7j)) = sign (ho(jj)) for all j = 1, 2, . . . , k — 1. By Theorem 14.31 the following inequalities hold 

< (i\ < /3i < 71 < (i\ < /3 2 < 72 < • • ■ < Mfe-i < Pk-i < 7fc-i < Mfc < /3fc, (5.16) 
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where j3j are the zeroes of the polynomial ho- From (I5.16[) it follows that ho does not change its sign on 
each interval (7j,/4 +1 ), j = 1, 2, . . . , k — 1, that is, sign (/lo(7j)) = sign (^o(^?+i)) f° r an J = 1, 2, . . . , fc — 1, 
in particular. Since h(/j,j+i) — and therefore fto(/i| +1 ) = — Mi^i(M?+i)) we have 

Poitf+i) = ' l o(Mi+i)5o(^ 2 +i) - Mi+ifro04fi)ffiOi+i) = ^o(Mj+i).9(-^+i)- 

The inequalities gr(— fJ.j+i) > 0, j = 1, . . . , fc — 1, proved above imply sign (i?o(M?+i)) = sign (^o(m?+i)) = 
sign (/io(7j)) = sign (po {lj ) ) • Consequently, the polynomial po has an even number of zeroes, counting 
multiplicities, on each interval {jj, /J/j-\-i), j = l,...,fc — 1. But the zeroes of the polynomials po and fti 
interlace, and /ii has no zeroes on the intervals (7^, /i^ +1 ), j = 1, . . . , k — 1, therefore po also has no zeroes 
on those intervals. So we proved the following inequalities 

^ < u t < 7j- < Mi+i for j = l,...,fc-l. (5.17) 

Let us consider the intervals (/j3,ujj), j = l,...,k and prove that the polynomial g(—y/u) has no 
zeroes on those interval. Suppose that it is not true. Then there exists a number /3 6 (fj3,0Jj) such that 
g(—VW) = 0. Then g (f3) = y/Pg\(fi) > 0, and from the equality 



Po 



08) = ho(p)go(0) + Phx(p)gi{p) = <?o(/3)Mv^). 



we obtain that sign(p (,3)) = sign(ft(y^)). It contradicts the equality sign(p (/3)) = — sign (h(\fj3)), 
which follows from (|5.17[) and from the equalities sign(po(0)) = sign(/io(0)) = sign(/i(0)) established 
above. 

Thus, we obtain that the polynomial g has real zeroes (en even number of zeroes, counting multiplicities) 
only on the intervals (—00, y^Wfe), (— p/., —y/uik-i), • • ■ , (— M2> — v 7 ^")' ( — Mi> 0)- Consequently, the following 
equalities hold sign (g(0)) — sign {g{—^j)) — sign (g l (— ->/wJ)) = sign(aj) > for all j = 1, 2, . . . , k. So we 

showed that if n = 2/ and p(0) ^ 0, then the function $(u) = — can be represented as in (I5.15[) with 

Po 
positive Uj and with exactly k positive poles, as required. 

Let now n = 2/ + 1 and p(0) 7^ 0. Then by definition, the polynomial p has only k zeroes in the closed 
right half-plane, all of which are positive: < fii < (12 < ■ ■ ■ < fJ-k- Moreover, p has an odd number 
of zeroes of p, counting multiplicities, on each interval (— /life,— /Ufe-i), ■■•, (— H3,— H2), (— P2,— Pi), and 
an even number of zeroes, counting multiplicities, on the intervals (— 00,— fXk) and (— /^i,0). All nonreal 
zeroes of p are located in the open left half-plane. 

As in the case n = 21, we represent the polynomial pas a product p(z) — g(z)h(z), where h is self- 
interlacing polynomial (of type I) of degree 2k— 1 whose positive zeroes are fi\, /12, • ■ • , fJ-k, and g is Hurwitz 
stable polynomial of degree 2(1 + 1 — k), which has an even number of zeroes, counting multiplicities, on 
each interval (—00, — Pfc), (—fik: ~ Pfc-i), ■ • •) ( — P2, — Mi)) ( — Hi, 0). In the same way as above one can prove 

that the function $ = — is an i?-function. The only difference between this case and the case of n = 21 
Po ^^^ 

is the number of positive poles of the function $. In fact, let us consider the function F defined in (J5.14I) . 

In this case, we have degpo = /if a\ ^ (or degpo = / — 1 if ai = 0), deg h\ = k — 1, deg g = / + 1 — k, 

therefore deg(/ii<?o) == /■ Thus, the number of poles of the function F is equal to (or greater than) the 

number of its zeroes. 

Let 

g{z) = 6 ^'- fe )+ 2 + Cl z 2 ^ +1 + . . . , 

h(z) = coz 2k - 1 + Cl z 2k - 2 + . . . 

Then the leading coefficient of the polynomial h\ equals Co, and the leading coefficient of the polynomial go 
equals 60 according to (|2.4j) — ()2.5[) . Consequently, the leading coefficient of the polynomial /iiffo is equal 
to co&o == o-o- Moreover, the leading coefficient of the polynomial po is equal to a\ oa 18 l ay, by (|2.5[) . so we 
have lim F(u) — or . The function F is an i?-function (that can be proved in the same way as 

u— s-oo ao ao 

18 If ai = 0. 
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above), so it is decreasing between its poles, and -F(O) = 7—7 > (see (14.201) ). Therefore, the minimal 

/ii(0)_ 
positive zero of the polynomial po is less than the minimal positive zero of the polynomial hi . 

If the condition (|5.3I) holds, then lim F(u) = < 0, so the maximal positive zero of po (the 

maximal positive zero of the function F) is greater than the maximal positive zero of the polynomial h\. 
This means that the polynomial po has exactly deg hi + 1 = k positive zeroes, so the function $ has exactly 
k positive poles. 

If the condition (|5.4[) holds, then either F(u) —t +0 as u — >• +oo if ai =0 or F(u) —} > as 

a ° 
u — > +oo if a\ ^ 0, so we have lim F(u) > 0. This implies that the maximal positive zero of hi is 

greater than the maximal positive zero of po, so the polynomial po has exactly k — 1 — deg hi positive 
zeroes. Consequently, $ has exactly k—1 positive poles. Thus, for n = 21 + 1 and p(0) ^ 0, we obtain that 

the function $ = — can be represented as in (|5.15l) with positive ctj and with exactly k (if (|5.3[) holds) or 

Po 
k — 1 (if (|5.4p holds) positive poles, as required. 

Let n — 21 and p(0) = 0. Then by definition, the polynomial p has exactly k zeroes in the closed right 
half-plane: = (Mi < jj,2 < . . . < (Ah- Moreover, p has an odd number of zeroes on each interval (—00, — Hk), 
(— /ifc, — /ife-i), . . . , (— /J>3, —(Jii), (— A*2, 0), counting multiplicities. All nonreal zeroes of p are located in the 
open left half-plane. 

As above, we represent the polynomial p as a product: p(z) — g(z)h(z) = po(z 2 ) + zpi(z 2 ), where 
h(z) = ho(z 2 )+zhi(z 2 ) is an almost self-interlacing polynomial 19 ! (of type I) of degree 2k whose nonnegative 
zeroes are /ii,/i2, • ■ • , /-tfc, and g(z) = go(z 2 ) + zgi(z 2 ) is a Hurwitz stable polynomial of degree 2(1 — k), 
which has an even number of zeroes, counting multiplicities, on each interval (— oo, —/ifc), (— fJ.k, — Mfe-i)) 
•■•, (-M3)-M2), (~M2,0). 

In the same way as above, one can establish that $ is an i?-function. However, po has a root at zero in 
this case. The polynomial ho also has a root at zero, and all its other roots are positive. Since deg h = 2k, 
we have deg/io = deg/ii + 1 = k. Besides, the zeroes of ho interlace the zeroes of hi, so all zeroes of the 
polynomial hi are positive. Consequently, the i?-function F defined in (|5.14|l has k—1 positive poles and 
I — k negative poles. It is easy to see that deg(go^i) =1 — 1 and degpo = I, so the maximal positive zero of 
the polynomial po (and of the function F) is greater than the maximal positive zero of the polynomial hi. 
Also the minimal positive of po is also greater than the minimal positive zero of hi, since F is decreasing 
between its poles and F(0) = 0. Thus, we obtain that po has exactly k—1 positive zeroes, since the zeroes 
of po interlace the zeroes of gohi. Consequently, the function $ has / — k negative poles, k—1 positive 
poles and a pole at zero. 

If n = 21 + 1 and the polynomial p(z) = po{z 2 ) + zpi(z 2 ) has a root at zero, then in the same as above, 
one can prove that $ is an i?-function with either exactly k—1 positive poles (if the condition (|5.3j) holds), 
or exactly k — 2 positive poles (if the condition (J5.4I) holds), since we have -F^O) =0 for the i?-function F 
defined in ([HTT3]) . □ 

Thcorem lS. 31 allow us to use properties of R- functions to obtain additional criteria for polynomials from 
the class of generalized the Hurwitz polynomials. In particular, we can generalize Hurwitz and Lienard 
and Chipart criteria. 

Let us consider the Hankel matrix S(&) = ||si+j||fj = o constructed with the coefficients of the se- 
ries (|2.13[) or (|2.18|) . From Theorems 15.31 11.18[ and 11.241 we obtain the following criteria for generalized 
Hurwitz polynomials. 

Theorem 5.4. The polynomial p defined in (|5.ip . with ai =/= if n = 21 + 1, is generalized Hurwitz if 
and only if 

r>i(*)>0, j = l,..., I. (5.18) 

Its order k equals: 

1) if n = 21, then 

k = l- V(l, 5i($), 5 2 ($), . . . , D P (*)), (5.19) 



19 See Section 1431 
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2) ifn = 21 + 1, then 

k = l- V(-«_i, 1, Si(*), S 2 (*)> . . . , D r ($)) + 1, (5.20) 

where r = I for 5/($) ^ 0, and r = I - 1 /o@ A($) = 0. Here s_i = — . 

XTie number of sign changes in the sequence 1, Z?i(<i>), Z?2 (<&),..., Z3 r (<I>) mwsi 6e calculated according 
to the Frobenius rule provided by Theorem \l.b\ 

Proof. According to Theorem 15.31 the polynomial p, with a\ ^ if n — 21 + 1, is generalized Hurwitz if 

r ji i 

and only if its associated function $ is an iZ-function with exactly I = — poles. But by Theorem II. 181 

this is equivalent to the inequalities (|5.18[) . 

If n = 21, or if n = 21 + 1 with (I5.3p . then by Theorem 15.31 the number of nonnegative poles of the 
function $ is equal to order of the polynomial p. From Theorem 1 1 . 241 we obtain (I5.19[) . which is equivalent 
to (J5~2TJ| if n = 21 + 1 and (JO) holds. 

If n. = 21 + 1 and aoai < 0, then by Theorem 15.31 the number of nonnegative poles of $ is equal to 
order of p minus one. Thus, by Theorem 11.241 we have 

k - 1 = i - v(i,5i(*), 5 2 ($), . . . , D r (*)), 

that is equivalent to (|5.20|) in this case. □ 

In the same way as above one can establish the corresponding theorem for polynomials of odd degree 
with a\ = 0. 

Theorem 5.5. Let the polynomial p of degree n = 2l-\-\ be defined in (|5.1[) and let a\ = 0. The polynomial 
p is generalized Hurwitz if and only if aoa^ < and 

£>i($)>0, j = l,...,l-l. (5.21) 

Its order k equals: 

k = l-V(l,D 1 ($),D 2 ($),...,D r ($)), (5.22) 

where r = I - 1 for A_i($) ^ 0, and r = 1-2 /o@ A_i($) = 0. 

The number of sign changes in the sequence 1, Z?i(<i>), D2 (<&),..., -D r (<I>) mwsi 6e calculated according 
to the Frobenius rule provided by Theorem \1.6\ 

Proof. In fact, let the polynomial p is generalized Hurwitz of order k. Then by Theorem l5.3l the function 

$ = — is an _R-function. In particular, this means that |degpo — degpi| ^ 1 (see e.g. [16]), therefore, 

Po 
degpi = I and degpo = Z — 1, since degp = 21 + 1 and a\ = by assumption. So by Theorems II. 181 and 15. 31 

the function $ has the form 

Hu) = ^ u+ a ^- a ^ + y^ 1 _ + y^^ (5>23) 

15 d j=i ■> j=k 3 

where — < 0, olj > for j = 1, . . . , I — 1, uij ^ 0, j = 1 . . . , k — 1, and i/j > for j = k, . . . ,1 — 1. The 

03 
inequalities f|5.21[) and the formula (J5.22I) now follow from Theorems 11.181 and 11.241 

Conversely, from the conditions (|5.21l) - (|5.22[) . a\ = and from the inequality aoa^ < it follows that 

the function $ has the form (|5.23[) by Theorems 11.181 and 11.241 Now Theorem 15.31 implies that p is a 

generalized Hurwitz polynomial of order k. □ 

From Theorems 15.41 and 15.51 and from the formulas (|2.14|) - (|2.15p and (|2.21|) we obtain the following 
theorem, which is an analogue of Hurwitz stability criterion for generalized Hurwitz polynomials. 

20 By Corollary 1131 if A(<t>) = 0, then i5;„ 1 (<S>) ^ 0. 
21 By Corollary H31 if D^^*) = 0, then S ; _ 2 (*) ¥" 0- 
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Theorem 5.6 (Generalized Hurwitz theorem). The polynomial p given in (|5.1[) is generalized Hurwitz if 
and only if 

A„_ x (p) > 0, A„_ 3 (p) > 0, A n _ 5 (p) > 0, . . . (5.24) 

The order k of the polynomial p equals 

fe = V(A„(p),A n _ a (p),...,l) if p(0)?0, (5.25) 

or 

fc = V(A B _a(p),A n _ 4 (p),.-.,l) + l «/ p(0) = 0. (5.26) 

//ere i/ie number of sign changes in (|5.25p and (|5.26p must 6e calculated according to the Frobenius rule 
provided by Theorem ] l.b\ 

Proof. Let the real polynomial p be of degree n such that the coefficient a\ is nonzero for n = 21 + 1. 
By Theorem l5.4[ the polynomial p is generalized Hurwitz if and only if the inequalities (j5.18[) hold. These 
inequalities are equivalent to the inequalities ()5.24|) according to the formulae (|2.14[) - (|2.15|l . 

If the polynomial p is of odd degree and ai = 0, then by Theorem l5.51 the polynomial p is generalized 
Hurwitz if and only if the inequalities (|5.2ip and — aoa3 > hold. Those inequalities are equivalent to the 
inequalities (|5.24j) according to the formulas (I2.21l) - (|2.22p . 

We now prove that order of the polynomial can be calculated by the formulae (|5.25p ~ (|5.26[) . 

Let n = 21, then from (|5.19[) and (|2.14p we obtain 

k = l- V(l, Di($), D 2 (*), . . . , fi r (*)) = I - V(l, -A 2 (p), A 4 (p), . . . , (-l) r A 2r (p)) = 

= V(l, A 2 (p), A 4 (p), . . . , A 2r (p)) + / - r. 

This is exactly (I5.25P for r = Z and exactly (|5.26p for r = I — 1. 

Let now n = 21 + 1 and a a! > 0. The formulas (|5.20p and (|2.14|l imply 

fc = j - v(-*_i, i, Si(*), S 2 ($), . . . , 5 r ($)) + i = 

= 1 + 1- V(-s_i, 1) + V(l, Di($), B 2 ($), . . . , D r ($)) = 

= I - V(l, -a^ 3 A 3 (p), a- 5 A 5 (p), . . . , (-l)- a r 2r - 1 A 2r+ i(p)) = 

= V(l,Ai(p),A3(p),A 6 (p),...,A 2r+ i(p)) + J-r, 

since a\ = Ai(p) > 0. One can see that (|5.27p implies (J5.25I) for r = I and (|5.26p for r = I — 1. 
Let n = 2l + l and a ai < 0. The formulae (j5.20|) and (|2.14p imply 

fc = i - v(-s_i, i, Bi(*), 5 2 ($), . . . , d p (*)) + 1 = 

= i-V(-3_ 1 ,l)+V(l,5 1 ($),B 2 ($),..., J D r ($)) + l = 

= 1 + / - V(l, -ar 3 A 3 (p), ar 5 A 5 (p), . . . , (-l) r a^ r - l A 2r+1 (p)) = 

= 1 + V(oi, a- 2 A 3 (p), a" 4 A 5 (p), . . . , ai " 2r A 2r+1 (p)) + 1 - r = 

= V(l, Ai(p)) + V( ai , ar 2 A 3 (p), ar 4 A 5 (p), • ■ ■ , ar 2r A 2r+1 (p)) + / - r = 

= V(l, Ai(p), A 3 (p), A 5 (p), • ■ • , A 2r+1 (p)) + / - r, 

since ai = Ai(p) < 0. Now from (f5T^5|) we obtain (fST2"5|) for r = I and (j5~2"6")l for r = I - 1. 



(5.27) 



(5.28) 
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Finally, let n = 21 + 1 and a x = Ai(p) = 0. The formulae (15T2"2"|) and ([^T|) - (l^^l) yield 

k = l-V(l,D 1 ($),D 2 ($),...,D r ($)) = 

Ab(p) A r (p) , lXr A 2r+3 ( P ) 



t-l-r + V(l J > -a ol) + v(l, 



A 3 (p)'A 3 (f))'---' v y A 3 (p) 

A 5 (p) A 7 (p) A 2r+3 (p)\ (5.29) 



A 3 (p)'A 3 (p)'---' A 3 (p) 
= / - 1 - r + V(l, A 1 {p), A 3 (p)) + V (A 3 (p), A 5 (p), A r (p), . . . , A 2r+3 (p)) = 

= i-l-r + V(l ) Ai(p)A3(p),A 6 ( P ),A 7 ( P ),...,A 2l . + 3( P )), 

where r = / — 1 for p(0) ^ 0, and r = I — 2 for p(0) = 0. Now (|5.29|) is equivalent to (|5.25[) for r = I — 1 
and to (pT2"ol) for r = 1-2. U 



Remark 5.7. Using the formulas (|2.14|) — (|2.15|) and (|2.21|) (|2.22[) . one can easily reformulate Theorem 15.61 
in terms of the minors rjj (p) defined in Definition 12.21 

Now from Theorems 11.281 and 15.31 we obtain the following theorem, which is a generalization of the 
famous Lienard and Chipart criterion of stability for the class of generalized Hurwitz polynomials. 

Theorem 5.8 (Generalized Lienard and Chipart theorem). The polynomial p given in (|5.1[) is generalized 
Hurwitz if and only if the inequalities (|5.24|) hold. The order k of the polynomial p equals 

k = v(a n ,a n -2,---,i)=v(a n ,a n -i,a n -3,...,l) if a n ^ 0, (5.30) 

oB 

k = v(a n - 2 ,a n -4, ...,1) + 1 = v(a n -i,a n -a,...,l) + 1> if a n = 0, (5.31) 

where v is the number of strong sign changes (see Definition \1.27\ ). 

Proof. The fact that the polynomial p is generalized Hurwitz if and only if the inequalities (j5.24[) hold 
follow from Theorem 15.61 Now we prove the formulae (|5.30[) - (|5.3ip . 

Let n = 21. By Theorem 15.31 p is a generalized Hurwitz polynomial if and only if the function 

Pi(u) _ aiu 1 ^ 1 + a 3 w'~ 2 + . . . + a 2 ;_ 3 u + a 2 i~i 



$(u) = 



Po(u) a u l +a 2 u l 1 + . . . + a 2 i- 2 u + a 2 i 



is an i?-function. The number of its nonnegative poles is equal to order k of the polynomial p. Let $ has 
exactly m positive polea 23 !. We prove that 

m = v(a n ,a„- 2 ,. ■■,!)■ (5.32) 

It will be equivalent to to the first equalities in (|5.30|) - (|5.31|) . because v is the number of strong sign 
changes, so v(a n , a n _ 2 , . . . , 1) = u(a n _ 2 , a„_4, . . . , 1) if a n = 0. Now Theorem II .281 implies 

m = v(a 2 i,a 2 i-2,- ■ ■ ,a ) = w(a 2 ;,a 2 ;_ 2 , . . . ,a ) + v(a a , 1) = w(a 2 ;,a 2 ;_ 2 , . . . ,a , 1), 

which is exactly (|5.32p . 

Let us now consider the function 

1 Pa(u) a u l + a 2 u l ~ l + . . . + a 2i - 2 u + a 2i 
y/ ( ii ) = = = . 

$(u) Pi{u) aiu 1 - 1 +a 3 u l ~ 2 + . . . + a 2l _ 3 u + a 2 i-i 

This function is also an i?-function by Theorem II .231 since $ is an i?-function with exactly k nonnegative 
poles, say ^ u\ < . . . , Wk- Note that $(u) > for u £ (u>k, +oo), so it can have at least fc — 1 and at 
most k positive zeroes. There can be only four possibilities: 



22 lia n = 0, then always a„—i ^ 0, since the polynomials pi andpo, the numerator and the denominator of the R- function <E>, 
have no common roots. 

23 Clearly, m = k if a n j^ 0, and m = k — 1 if a n = 0. 
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I. 3>(0) = < 0. Then v(a2i,a2i-i) = 1, and $ has no zeroes on the interval (0,oji), since it is 

a 2 i 
decreasing on this interval and $(w) — > — oo as u /"■ uj\. Therefore, $ has exactly k — 1 positive 

zeroes, so $ has exactly fc — 1 positive poles. By Theorem 11.281 and Remark 11.291 applied to the 

function *, we have 

fc- 1 = v(a 2 i-i,a 2 i- 3 ,...,ai) = -1 + v(a 2 i,a 2 i-i) + v(a 2 i-i,a 2 is,--. ,ai,l) + u(ai,l) = 

= -1 + U(02J, 02J-1, 02J-3, • • • , Oil 1)) 

which is exactly 

fc — w(a„,a„_i,a„_3, . . . , 1). (5.33) 

Here we used the inequality Ai(p) = ai > 0, which follows from the inequalities (|5.24[) for n = 21. 

II. <1>(0) = — — > 0. Then v(a 2 i,a 2 i~i) — 0, and $ has exactly fc positive zeroes, so \I/ has exactly fc 
a 2 i ^^ ^^ 

positive poles. Thus, from Theorem II .281 and Remark 11.291 applied to the function "J, we obtain 

k = v(a 2 i-i,a 2 i-3, . . . ,Oi) = w(a 2 /,a 2 ;-i) + v(a 2 i-i,a 2 i-3, • ••>ai) + w(ai, 1) = 

= U(02i,a2j-l,02i-3i • • • i aii 1), 

which is exactly (|5.33p . 

III. 3>(0) = 0. Then \& has a pole at zero, that is, a 2 i-i = 0. Since ^> is an i?-function, it has positive 

residues at each its pole, so lim u^'(u) = > 0. Thus, v(a 2 i,a 2 i-3) — v(a 2 i,a 2 i-i,a 2 i-3) = 1. 

«-i-o a 2 i-3 

As in Case I, $ has exactly fc — 1 positive zeroes, so the function \1/ has exactly fc — 1 positive poles. 
Therefore, by Theorem 11.281 and Remark 11.291 applied to the function 'J, we have 

fc- 1 = v(a 2 i-3,a 2 i- 5 ,...,ai) = -1 +v(a 2 i,a 2 i-i,a 2 i- 3 ) +v(a 2 i-3,a 2 i- 5 ,...,ai) +u(oi,l) = 

= -1 + w(a 2 ;,a 2 ;-i,a 2 ;-3, ■ • ■ ,ai, 1), 
which is precisely (|5.33p . 

IV. \P(0) = 0. Then p(0) = a n — a 2 i = 0, and the function $ has a pole at zero, so it has exactly fc — 1 
positive zeroes. Thus, the function * has exactly fc — 1 positive poles, and by Theorem 11.281 and 
Remark 11.291 applied to the function ^, we have 

fc - 1 = v(a 2 i-i,a 2 i-3,- . . ,Ol) = w(a 2 ;_i,a 2 ;_3,. ..,Ol,l), 

therefore, 

fc = u(o„_i,a„_3, ...,1) + 1. (5.34) 



Let now n = 21 + 1 . By Theorem 15.31 p is a generalized Hurwitz polynomial if and only if $ is an 
i?-function with exactly I zeroes and I or I — 1 poles. If fc is order of the polynomial p, then $ has exactly 
fc nonnegative poles for a± > and exactly fc — 1 nonnegative poles for ai sj 0. 

At first, we establish the formula (|5.32[) . where m is the number of positive poles of the function $. As 
we mentioned above, this formula is equivalent to the first equalities in (|5.30p - (|5.3ip . 

If ai > 0, then from Theorem 11.281 applied to the function $ we obtain 

m = v(a 2 i + i,a 2 i-i, . . . ,oi) = v(a 2 i + i,a 2 i~i, . . . ,ai, 1). 

This is exactly (jQ2"]) . 

Let ai ^ 0. As we established above, if a\ — 0, then 03 < 0. Therefore, for a\ ^ 0, we always have 
v(ci3,ai) = ^(03, 01, 1) — 1. Thus, Theorem II .281 implies 

m = v(a 2 i + i,a 2 [-i, . . . ,Oi) = w(a 2;+ i,a 2; _i, . . . ,ai, 1) - 1, 
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which is again (I5.32J) . since for ax $J 0, m = k — 1 if p(0) ^ 0, and m = fe — 2 if p(0) = 0. 

As above, consider the function 

_ 1 p {u) _ axu 1 + a^u 1 " 1 + . . . + a 2 i-iu + q 2 ; + i 

$(u) pi (u) a u z + a 2 u /_1 + . . . + CL21-2U + a 2 i 

This function is also an R- function by Theorem ll.231 As above, there can be only four possibilities: 

I. $(0) = < 0, so v(a2i+i,a,2i) = 1. If a\ > 0, then by Theorem 15.31 $ has exactly k positive 

a-21+i 

poles, say < u>\ < ... < w*. Since $(u) — > — > as u — > +00 and $ is decreasing between its 

a\ 
poles, the function $ has no zeroes on the interval (w^, +00). Also $ has no zeroes on the interval 

(0,Wi), since $(0) < by assumption. Thus, $ has exactly fe — 1 positive zeroes, so $ has exactly 

fc — 1 positive poles. 

If a\ < 0, then by Theorem 15.31 the function $ has exactly k — 1 positive poles, say < wi < 
. . . < Wfc_i. Since $(w) is decreasing between its poles, it is negative for sufficiently large positive u. 
Moreover, <b{u) /* +00 as u \ OJk-i- Thus, $ has a zero on the interval (uik-i, +00), but it has no 
zeroes (0, Wi), since $(0) < by assumption. Consequently, $ has exactly k — 1 positive zeroes, so 
\1/ has exactly k — 1 positive poles. 

Thus, we showed that the function ^ has exactly k — 1 positive poles regardless of the sign of 01. 
Now from Theorem II .281 and Remark 11.291 applied to the function Vt, we obtain 

k— 1 = v(a 2 i,a 2 i-2,- . .,00) = — 1 +^(021+1,021) + v(a 2 i,a 2 i-2, ■ . .,00, 1) = 

= -1 + V(02J+1, 02J, 02J-2, • • • , Oo, 1), 

which is exactly (I5.33[) . 

Ct97 

II. $(0) = > 0, so v(a2i+i, a>2i) — 0- I n the same way as above, one can show that the function \1/ 

a-21+i 
has exactly k positive poles regardless of the sign of a\. Thus, from Theorem 11.281 and Remark ll.291 

applied to the function ^, we obtain 

k = v(a 2 i,a 2 i-2, ■ ■ ■ ,Oo) = u(o2/+i,02/) + v(a 2 i,a 2 i-2, ■ ■ ■ ,Oo, 1) = 

= V(02l+1, 021,021-2,. ■ ■ ,O0, 1). 

So the formula (|5.33j) is also valid in this case. 

III. $(0) = 0. Then ^ has a pole at zero, that is, 021 = 0. Since "J is an i?-function, it has positive residues 

at each its pole, so lim u^(u) = > 0. Thus, v{a,2i+i, 0,21-2) = v(o2i+i, 021, 0,21-2) = 1. As 

u^o a 2 /-2 

above, one can show that VP has exactly fe — 1 positive poles regardless of the sign of ai. Therefore, 
by Theorem 11.281 and Remark 11.291 applied to the function "J, we have 

k - 1 = v(a2i-2,a.2i-4,- ■ ■ ,a ) = -1 +V (021+1, 0,21, 021-2) + v(a 2 i-2, a 2 i-A, ■ ■ ■ ,a , 1) = 

= -1 + v(a2i+i, 021, 0,21-2, ■■■,ao, 1), 
which is precisely (J5.33I) . 

IV. ^(0) = 0. Then p(0) = a n = a 2 ;+i = 0, and the function $ has a pole at zero. In the same way as 
above, one can prove that the function \1/ has exactly fe— 1 positive poles regardless of the sign of 01, 
so Theorem 11.281 and Remark 11.291 applied to the function "J imply 

k - 1 = v(a 2 i,a 2 i-2, ■ ■ ■ ,00), 
which is equivalent to (|5.34p . since ao > 0. 
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D 



This theorem implies a necessary condition for polynomials to be generalized Hurwitz. This condition 
plays a role of Stodola's necessary condition for Hurwitz stable polynomials (Theorem I3.2[) . 



Corollary 5.9. If the polynomial 

p(z) = a^z 71 + a\Z n ~ + • • • + a n , oi, . . . , a n € R, oq > 0, (5.35) 

is generalized Hurwitz (of type I) of order k, then 

fc = «(a„,a„_2,...,l) = u(a„,a„_ 1 ,a n _ 3 , ...,1) if p(0) ^ 0, 

or 

k- 1 = v(a„, a„_ 2 , • • • , 1) =t>(a„,a„_i,a„_3, •••,!) «/ p(0) = 0, 

This corollary implies Theorem 13.21 (Stodola's theorem) for /c = 0, and it implies Theorem 14.101 for 

' n + ll 
fc = — - — and a n 7^ 0. 

Note that G. F. Korsakov [21] made an attempt to prove Theorem 15.81 However, his methods did 
not allow him to prove the whole theorem. He proved only sufficiency, that is, he proved that if the 
inequalities (|5.24[) hold for the polynomial (|5.35[) . then p is generalized Hurwitz, and the number k of its 
positive zeroes can be found by the formulae 

• for n = 21, 

fc = w(a„,a„_ 2 , ■ • • ,a 2 ,a ); (5.36) 

• for n = 21 + 1, 

k = v(a n , a„_ 2 , • • • , a 3 , oi, a ). (5.37) 

Since in [21] . it was assumed that a n ^ 0, the number k in (I5.36[) - (|5.37[) is exactly order of the polynomial p. 
Nevertheless, the generalized Hurwitz polynomials with a root at zero are also partially described in [2Tj . 
In fact, there was proved that if the following inequalities hold for the polynomial p defined in (j5.35[) 

A n (p) > 0, A„_ 2 (p) > 0, A n _ 4 b) > 0, . . . , (5.38) 

then j? has exactly k (defined by (jn3S])~(IS3Zl)) zeroes on the positive real half-axis, and other its zeroes are 
located in the open left half-plane. Generally speaking, the polynomial p satisfied the inequalities (15.38)) is 
not a generalized Hurwitz polynomial. However, if p satisfies (I5.38|) . then the polynomial q(z) = zp(z) of 
degree n + 1 is a generalized Hurwitz polynomial of order fc + 1. Indeed, it is easy to see that the following 
equalities are true: 

Ai(q) = Ai{p), i = l,...,n. 

So the inequalities (I5.38|) for the polynomial p of degree n are exactly the inequalities (|5.24[) hold for the 
polynomial q of degree n + 1, so q is a generalized Hurwitz polynomial having k + 1 nonnegative zeroes. 

At last, we establish a relation between generalized Hurwitz polynomials and continued fractions of 
Sticltjes type (if any). 

Theorem 5.10. // the polynomial p of degree n defined in (J5.35J) is generalized Hurwitz (of type I) of 
order k and A n _2j(p) 7^ 0, i = 1, . . . , I, where Ao(p) = 1, then the function $ has the following Stieltjes 
type continued fraction expansion: 

<, {u) = eM = C0 + 1 , T = { <* vmfo, (5 . 39) 

V ^ Po(u) 1 1 c 2l ^u if p(0 =0, y ' 

C\U -I '■ 



C 2 + 

C3U 



1 
f 
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where Cq = if n is even, and Cq 7^ if n is odd, and 

c 2j -i>0, j = l,...,l. (5.40) 

Moreover, if p(0) 7^ 0, then the number of negative coefficients C2j, j = 0, 1, . . . , I, equals k, and if p(0) = 0, 
then the number of negative coefficients cij, j — 0,1, ... ,1 — 1, equals k — 1. 

Conversely, if for given polynomial p of degree n, its associated function $ has a Stieltjes continued 
fraction (|5.39[) f)5.40[) . thenp is a generalized Hurwitz polynomial (with A n -2j(p) =/= 0, i = 1, . . . , I), where I 
defined in (|2.2j) . Its order equals the number of negative coefficients C2j , j = 0, 1, . . . , /, if p(0) 7^ 0, or the 
number of negative coefficients C2j, j = 0, 1, . . . , I — 1, plus one if p(0) = 0. 

Proof. The theorem follows from Theorem 15.61 and from the formula? p,13[) - (l3.14p □ 

Theorem 15.101 does not cover the case of ai =0. The following theorem fills this gap. 

Theorem 5.11. // the polynomial p of degree n = 21 + 1 defined in (|5.35[) with a± — is generalized 
Hurwitz (of type I) of order k and A2j+i(p) =/= 0, i = 1, ...,/ — 1, then the function $ has the following 
Stieltjes continued fraction expansion: 

®W = —— = -c- lU + c Q + , T=< (5.41) 

Po{u) 1 C21-3U if p(0) =0, 

c\u H v 

c 2 + : 



c 3 u 



1 

T 



where cq € R, and 



c 2l -i>0, j = 0,l,...,l-l. (5.42) 

Moreover, if p(0) 7^ 0, then the number of negative coefficients C23, j = 1, ...,/ — 1, equals k — 1. But if 
p(0) = 0, then the number of negative coefficients c%j, j = 0, 1, . . . , I — 2, equals k — 2. 

Conversely, if for given polynomial p of degree n = 21 + 1, its associated function $ has a Stieltjes 
continued fraction (|5.41j) - (|5.42j) . then p is a generalized Hurwitz polynomial (with A2j+i(p) 7^ 0, i = 
1, . . . , I — 1 and Ai(p) = a\ =0). Its order equals the number of negative coefficients C2j, j = 1, ■ ■ ■ ,1 — 1, 
plus one if p(0) 7^ 0, or the number of negative coefficients C2j, j — 1, . . . ,1 — 2, plus two if p(0) = 0. 

Proof. Let the polynomial p of degree n — 21 + 1 with a\ — be generalized Hurwitz (of type I) of order k 
and A2j+i(p) 7^ 0, i = 1, . . .,1 — 1. By Theorem 15.51 and by formula? (J2.21I) . we obtain 0003 < and 

r d,-(*)>o, j=i,...,i-i, 

\ £>j($)^0 j = l,...,l-2. 
It is easy to see that for the function F(u) :— $(w) + c_iit, where c_i = , the following equalities hold 

{ D J F ^ D J^ J = l,...,l-1, (5.44) 

\ D j (F)=D j ($)^0, ' 

so the function F has a Stiltjes type continued fraction (jl.331) with C2j-i > 0, j = 1, . . . , I — 1, that follows 
from (|5T4l?|) - ([5T34l and (fOSl) . Also from (J5T2"2"1) . (|5^3l) - (|5Tl4l and pTBTl) we obtain that the number of 
negative coefficients C2j, j = 1, ...,/ — 1, equals fc — 1 if p(0) 7^ 0, or k — 2 if p(0) = 0. 

The converse assertion of the theorem follows from the formulae (|5.43p - (l5.44[) . from Theorem 11.361 

applied to the function F(u) := $(u) u, and from Theorem 15.51 □ 

03 



Theorem 11.221 implies the following theorem. 
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Theorem 5.12. Let the polynomial p of degree n ^ 2 as in (|5.35p be generalized Hurwitz polynomial. 
Then all polynomials 



n-2j 



are also generalized Hurwitz. 



n — i 


( 


n — i 


2 


2 



+ j-l)a i z n - 2j - i , k = l, 



5.2 Application to bifurcation theory. Polynomials dependent on parameters. 

Let us consider a system of ordinary differential equations dependent on a real parameter, say a: 



x(t) = F(x(t), a) = A(a)x(t) + o(x(t)), 
x(t ) = x 



F(x), x e 



(5.45) 



where the n x n matrix A(a) is the linear part of the function F(x,a), and F(0,a) = 0. Suppose that 
for some value of the parameter ao all the eigenvalues of the matrix A(ao) are located in the open left 
half-plane of the complex plane. In this case, the stability theory says [20, 2HJ [TB] that the zero solution 
of the system (J5.45I) is (asymptotically) stable, that is, all the solutions of the system (J5.45I) . whose initial 
conditions xo are sufficiently small, tend to zero as t — > +oo. 

When the parameter a arrives at its critical value, say a\ the following two basic cases can occur, 
generically: either a pair of complex conjugate eigenvalues (of multiplicity one) of the matrix A(a{) 
appear on the imaginary axis or an eigenvalue (of multiplicity one) becomes 0. In the former case, for 
a = a* + e, where e > is sufficiently small, there appears a small-amplitude limit cycle branching from 
the zero solution (a fixed point) of the system (|5.45p . but the zero solution loses its stabilit}0. This type 
of bifurcation is called Hopf bifurcation [251 HH HE] • in hydrodynamics [42 , such type of bifurcation is 
called an oscillatory loss of stability or oscillatory instability. In the case, when the matrix A(a\) has 
all the eigenvalues in the open left half-plane except one, which is equal to zero and is of multiplicity 
one, new fixed points of the system (|5.45j) branch from the zero fixed point. In hydrodynamics [42], such 
type of bifurcation is called a monotone loss of stability or monotone instability. To determine the type of 
bifurcation of a fixed point of a given system is a very important problem in hydrodynamics and mechanics 
(see [551 S2] and references there) . 

Let us consider the characteristic polynomial p(z, a) = det(zE — A(a)) of the matrix A(a). By 
assumption p(z,cto) is Hurwitz stable. Recall that a\ > a® is the critical value of the parameter a such 
that the zero solution of the system (|5.45j) loses its stability. The monotonicity principle takes place if and 
only if the polynomial p(z, a*) is a generalized Hurwitz polynomial, that is, the inequalities (J5.24J) hold. 
This fact follows from Theorem 15.01 Moreover, the polynomial p(z,a) is generalized Hurwitz of order 
(Hurwitz stable) for a < ai*, and p{z, a*) is a generalized Hurwitz polynomial of order 1. 

Now we describe further changes of zero location of the polynomial p{z, a) when a goes beyond a\ 

\ n + 1" 
until p(z, a) becomes self-interlacing. Let {a^, a^, . . . , «„}, H™^ — z — 

of the parameter a at each of which order of the generalized Hurwitz polynomial p(z,a) increases. This 
means that p(z, a*) has a simple root at zero, and for a > a* sufficiently close to a*, this root becomes a 
positive simple root of p{z,a) with the minimal absolute value (among real roots of p(z,a)). Therefore, 
when a transfers ot* m , the polynomial p(z, a) has exactly m positive simple zeroes for a > a* m sufficiently 
close to a* m . Denote the minimal positive zero of p{z, a) by /i m . Obviously, if the difference a — a* m > 
is sufficiently small, then the number fi m > is close to zero, so p(z, a) has no zeroes on the interval 
(— /j, m , 0). Let a increase. One more simple positive zero of the polynomial p(z, a) can appeaa 25 !. when two 
zeroes of a pair of complex conjugate zeroes of p(z, a) (all its nonreal zeroes lie in the open left half-plane) 
face each other on the interval (— fj, m , 0) to become a real zero of p(z, a) of multiplicity two. Further, when 
a grows, this zero splits into to real simple roots on the interval (— ju m , 0). One of those roots goes toward 
zero and crosses it, when a transfers through the next critical value a* n+1 of the parameter a. Thus, for 



be the next critical values 



4 Here we assume that for a > aj, the matrix A(a) has an eigenvalue with a positive real part. 

\n + 1" 
3 It is possible only if m < 
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sufficiently small difference a — a* n+1 > 0, this zero becomes the smallest positive zero p m +i- The second 
zero of the mentioned pair of zeroes remains on the interval (— p m , —fJ-m+i)-, since it cannot cross the value 
— fi m and cannot go away from the real axis. 

Let us show that there are no any other scenarios for positive zeroes to appear if p(z, a) is generalized 
Hurwitz for all a. In fact, suppose that for a > a* m , some negative zeroes of p(z, a) come into the interval 
(—fj, m , 0) from the real axis, that is, from the interval (— oo, —fi m )- But in this case, the polynomial p(z, a) 
will have zeroes ±fJ, m at some value of a, so by the Orlando's formula [12j US], we have A„+i(p) = 0. 
According to Theorem 15. 6[ this contradicts the assumption that p{z, a) is generalized Hurwitz. 

Thus, when a increases from ao, at each interval (— pj, — Pj+x) the polynomial p(z,a) has an odd 

number of zeroes. One of those zeroes appears on (— pj, — Pj+i) as it was described above, while the 

other zeroes appear from complex conjugate pairs. The maximal number of positive zeroes of p(z,a) 

\n + ll 
equals . In this case, p{z, a) becomes self-interlacing. 

So the monotonicity principle takes place for the system (|5.45[) if and only if for the first critical value a\ 
of the parameter a the characteristic polynomial p(z, a\) — det(zE — A(a\)) is generalized Hurwitz, that 
is, the inequalities (I5.24p hold. 
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